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Abstract. The control of internal and motional quantum degrees of freedom of laser cooled 
trapped ions has been subject to intense theoretical and experimental research for about three 
decades. In the realm of quantum information science the abihty to deterministically prepai'e 
and measure quantum states of ti'apped ions is unprecedented. This expertise may be employed 
to investigate physical models conceived to describe systems that are not directly accessible 
for experimental investigations. Here, we give an overview of current theoretical proposals 
and experiments for such quantum simulations with trapped ions. This includes vaiious spin 
models (e.g., the quantum transverse Ising model, or a neural network), the Bose-Hubbard 
Hamiltonian, the Frenkel-Kontorova model, and quantum fields and relativistic effects. 
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1. Introduction 

Atomic ions confined in an electrodynamic trap provide us with individual quantum systems 
whose internal and external degrees of freedom can be controlled to a considerable degree by 
the experimenter. A strong motivation for early experiments with individual trapped ions arose 
from their potential use as frequency standards (e.g., [Dehmelt 1981] and references therein). 
Today, by using trapped ions as frequency standards unsurpassed accuracy and precision has 
been reached (e.g., [Rosenband et al. 2008] and references therein), and in addition, trapped 
ions are used for a wide range of investigations into fundamental questions of physics, for 
instance in the field of quantum information science. This special issue gives an account of 
the many facets of research with trapped ions. 

Controlling motional and internal states of trapped ions by letting them interact with 
electromagnetic fields, ranging from uv Ught to rf radiation, has lead to the reaUzation of 
experiments that previously (reaching well into the second half of the 20th century) were 
only conceivable, by most physicists, as Gedankenexperiments. A small (and unsystematic) 
sample from such fascinating experiments include the trapping and visualization of a single 
ion [Neuhauser et al. 1980], the laser cooling of a single ion to its motional ground state 
[Diedrich et al. 1989], the deterministic generation and analysis of Schrodinger cat and other 
entangled states of massive particles, [Monroe, Meekhof, King & Wineland 1996, Leibfried 
et al. 2005, Haffner et al. 2005, Blatt & Wineland 2008], or the teleportation of ions [Riebe 
et al. 2004, Barrett et al. 2004]. 

While building a universal quantum computer still poses formidable experimental 
challenges, in particular because of the low error thresholds required for fault-tolerant 
computing it appears that useful quantum simulations are more amenable to experimental 
efforts. 

There are two common usages of the term "quantum simulation": (i) the simulation 
of static and dynamic properties of quantum systems (e.g. atoms, molecules, condensed 
matter systems) using classical computers, and (ii) the simulation of the properties of one 
quantum system by means of controlling and observing another quantum system. Here, we 
are concerned with the second type of research. 

The purpose of a quantum simulation is to use a well-understood physical A system to 
simulate the static properties and the dynamics of another system B that is difficult or, for all 
practical purposes, even impossible to investigate experimentally [Feynman 1982, Feynman 
1986, Jane et al. 2003]. In order to be able to simulate B by observing A, both physical 
systems need to be described by the same underlying mathematical model. 

In this article we give an overview of concrete proposals for performing experimental 
quantum simulations with laser cooled trapped ions, and of experimental work unmediately 
relevant for using trapped ions as a quantum simulator. 

2. Spin models 

Interacting spins play an important role in many physical models in a variety of research fields 
of physics, for instance, in models that are employed to describe phenomena in condensed 
matter physics (e.g., [Schollwock et al. 2004] and references therein). Also, in recent years 
the connection between entanglement on the one hand, and static and dynamic properties 
of many-body systems on the other hand, for instance, near quantum phase transitions has 
attracted much interest (e.g., [Osterloh et al. 2002, Amico et al. 2008, Plenio & Virmani 2007] 
and references therein). 
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Probably the simplest approach to modeUng interacting spins is the classical Ising model 
[Lenz 1920, Ising 1925] where the energy of a system of locaUzed particles is given by 

E=-JY, SiSj (1) 

and each particle (spin) maybe in two possible states indicated by Si = ±1. The summation 
includes the interaction between nearest neighbours only, indicated by {i > j). Beyond 
interacting spins, the Ising model is often used to approximate systems where the constituents 
may be in one of two possible states, describing, for example, the tunneling between two 
potential wells or the occupation of a lattice site with two types of particles [Brush 1967]. 
Ising already showed that a phase transition between an ordered state and a disordered state 
does not take place in this model in one dimension at non-zero temperature, whereas in 
two dimensions such a phase transition is possible as shown by Onsager [Onsager 1944]. 
However, extending the range of interaction beyond nearest neighbours gives rise to phase 
transitions even in one dimension [Gitterman & Halpern 2004]. 

A more sophisticated model of the interaction between locaUzed spins [Heisenberg 1928] 
that, in addition, are exposed to an external local magnetic field -Bj uses the Hamiltonian 

Hs = —— JijSi ■ Sj — Si - Bi . (2) 

i>j i 

Here we assume that all lattice sites are populated by particles with spin Si with identical 
gf-factor and eigenvalues of S^ that are multiples of U/2 (hb is the Bohr magneton). No 
assumption has yet been made on the physical nature and range of the spin-spin coupling 
which often is caused by a dipole-dipole or exchange interaction. Therefore, it is useful to 

(a) 

explicitly parametrize a possible anisotropy of this coupling by introducing the constants J-j ' , 
with a = x,y^z^o that 

^i"'=-^ E Y.j'S'st'sf^-Y.Si.B[, (3) 

a=x,y,z i>j i 

where the scaled local fields are defined as = ^^Bi. 

When describing physical phenomena that occur in condensed matter, for example, 
concerning magnetism, this Hamiltonian is often subject to idealizations with regard to its 
dimensionality, the range of interaction, its (an-)isotropy, the uniformity of couphng constants 
Jij, the magnitude and direction of additional (local or global) fields, and regarding the finite 
size of a particular system under consideration. 

The Hamiltonian of the quantum transverse Ising model is obtained from ifi"^ for spin- 
1/2 particles (i.e., S = {h/2)a) by setting J-^^ = Q = J^r'^ and then including further 

restrictions: Specifically, the local fields Bi are replaced by a global field B = {B, 0, 0)^ 
pointing in the x— direction, thus neglecting inhomogeneities in the physical system to be 
modeled [Igloi et al. 1993, Platini et al. 2007]. Additional approximations often made are 
that the interaction between spins is restricted to nearest neighbours, and uniform coupling 
jj^j^ = Jyi,j is assumed. This then leads to 

HTi = -P^ata]-B'^J2<- ^4) 

{i>3) i 

The quantum transverse Ising model is a prominent example for a spin model [Pfeuty 1970] 
that has been the starting point for many theoretical investigations (e.g., [Das & Chakrabarti 
2005] Part I and references therein). 
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2.1. Simulation of spin models with trapped ions 

Using two internal states of a trapped ion as an effective spin- 1/2 with the ions interacting 
via the Coulomb force allows for engineering spin couplings by letting the ions interact, 
in addition, with external dynamic and static fields. Initial theoretical and experimental 
work was mainly concerned with creating effective interactions between two spins to be 
used as 2-qubit quantum gates [Cirac & Zoller 1995, Monroe, Meekhof, King, Itano & 
Wineland 1995, Molmer & Sorensen 1999, Milburn 1999, Milburn et al. 2000, S0rensen & 
M0lmer 2000, Jonathan et al. 2000, Jonathan & Plenio 2001, Cirac & Zoller 2000, Sasura & 
Steane 2003, Schmidt-Kaler et al. 2003, Leibfried et al. 2003] 

2.1.1. Direct spin-spin interaction In this section we concentrate on couplings of the type 
given in equation |3] that, for example, allow to study a phase transition between ordered and 
disordered phases of an ion crystal as will be outlined in section |Z2l The question we want to 
consider first is. How can such a spin-spin coupling be obtained in an ion trap and how can it 
be engineered to simulate a specific Hamiltonian? 

If we consider a laser cooled Coulomb crystal of singly charged ions confined in a linear 
trap such that they arrange themselves in a 1-d string, then the relevant length scale for the 
inter-ion distance is 

C = (eyA^eomulf/^ , (5) 

where m is the mass of one singly charged ion, e the elementary charge, and vi is the 
angular vibrational frequency of the center-of-mass axial mode of the collection of ions 
[Steane 1997, James 1998]. We have C = 1.5 x l{r^{Mpy^/'^ m with M in atomic mass 
units and / = i^i/27r in units of MHz. Thus, C, typically amounts to a few micrometers. For a 
linear ion string the inter-ion spacing is 5z « C,2N^^-^^ [Schiffer 1993, Dubin 1993] which 
gives 3.7 /im for ^^^Yb+ and 9.9 /im for ^Be+ assuming two ions confined in a potential 
characterized by / = 1 MHz. 

In principle, a direct spin-spin interaction may take place via a magnetic moment 
associated with the (pseudo-)spin of the ions or via an exchange interaction. However, as will 
be recapitulated now, the inter-ion spacing renders this interaction negligible. It is nevertheless 
instructive to have a look at the dipole-dipole interaction, since it may help to gain an intuitive 
understanding also of the indirect spin-spin coupling treated in section |2".1.2| 

For the moment, for the sake of obtaining an intuitive picture, we consider two spins 
with their relative alignment fixed, and exposed to local magnetic fields such that the 

uncoupled Hamiltonian reads Hq — (?i/2)(ajicrp'' + lj2<j'"2^) (1^1.2 are the angular Larmor 
precession frequencies). Treating the magnetic dipole-dipole interaction as a perturbation to 
first order, its magnitude is determined by the expectation value 

(si,S2| Vmd |si,S2) = siS2fe{Sz) = ±?iAw (6) 

where si^2 = ±1 and 

TA Mo 71 72 2in\ (^) (^) 

Vmd - - Ti^^(3cos {9) ~ l)a\ >a\ 

= fe{Sz)<j['^4'^ . (7) 

Here, the gyromagnetic ratio of spins 1 and 2 is denoted by 71^2, Mo indicates the permeability 
of empty space, and 9 is the angle indicating the relative orientation of the two spins. Equation 
|6]tells us that the energy of the two-spin systems depends on the relative orientation between 
the two spins, that is, flipping either one of the two spins results in a change of energy. The 
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physical reason is that the energy of a given spin depends not only on the imposed local field 
B\ 2 , but also on the field generated by the second spin at the location of the first one, which 
in turn depends on its orientation (compare figure[T]and figure|2]l. 




Figure 1. An artist's illustration of the coupling between two spins via a dipole-dipole 
interaction. Two classical magnetic moments 7^5^ fixed in space are exposed to individual 
external magnetic fields Bi {i = 1,2) that determine their respective energy for a given 
orientation (parallel or anti-parallel to Bi). Spin 2 experiences in the presence of spin 1, 
in addition to B2, an additional field that depends on the relative orientation between the two 
spins. Thus the energy of spin 2 is shifted depending on the orientation of spin 1 (and vice 
versa). The spatial wavefunctions of the two spins are indicated by >I>i(ri). 

For two electron spins (71,2 ~ e/rrie with mass mg) and assuming 6* = 0, we find 
Aw « 27r X 2.6 x lO^'^^Sz^^ Hz with 5z measured in meters. Therefore, for a typical inter-ion 
separation, 5z — 5/im the direct magnetic spin-spin interaction amounts to Aw « 27r x 0.2 
mHz and is negligible in typical ion trap settings. This estimate is based on a magnetic 
dipole-dipole interaction of ions in their ground state. If bringing the ions to a highly excited 
Rydberg state, one would expect an appreciable interaction between their electric dipoles 
[Miiller et al. 2008, Urban et al. 2009, Gaetan et al. 2009]. 

In order for the exchange interaction to play a non-negligible role, an overlap of the wave- 
functions describing the external, motional degrees freedom of neighbouring ions would be 
required, or more precisely, the exchange integral / [$ri(n.)$n+i(?^2)] ['^'n (^2)^^+1 (n.)] 
would have to have a non-negligible magnitude. Here, the spatial wavefunction of ion number 
n (n + 1) located at position ri (^2) is <&n(n.) ($ri+i(?^2)) and the integration is carried out 
over all space. However, the root mean square extension of the ground state wavefunction of 
a trapped ion is given by Azi — yjh/2mvi ~ 7.1 x 10^^ x (A//)^^/^ m In the case of a 
^Be"*" (^^^Yb"*") ion, one obtains Azi — 24 nm (5.4 nm) for / = IMHz which typically is 
three orders of magnitude smaller than the inter-ion separation 5z. Therefore, the exchange 
interaction can be neglected, too. 



V^ {E, + E2) - flA(i) 
V2 (E, - E2) + ^Aco 



V2 i-E, + E^) - /zAco 
Yi i-E, - E2) + fiAd) 

r orientation. Tlie quantum numbers 
ncies of non-interacting spins are uii 

2.1.2. Spin-spin interaction mediated by the Coulomb force Thus, a direct spin-spin 
interaction between trapped ions will not be sufficiently strong for a useful implementation of 
the spin models described above. Instead, we now look for an interaction that is mediated by 
some other physical mechanism. 

Physical picture The applied trapping potential and the Coulomb forces between the ions 
determine their equilibrium positions. Now, by applying a suitable external field, an additional 
force is exerted on each ion whose magnitude and/or direction depends on the orientation of 
the spin of the ion and varies spatially (i.e., depends on the exact location of the ion). If 
this state-dependent force acts in addition to the trap and Coulomb forces, then a spin flip of 
ion 1 changes this force and, as a consequence, this ion changes its equilibrium position. Via 
Coulomb repulsion, this in turn leads to a shift of the equilibrium position of the neighbouring 
ion 2 that now finds itself exposed to a field of different magnitude and thus changes its 
internal energy. 

This additional field may be a static magnetic field with spatially varying magnitude. The 
direct interaction between magnetic dipoles (where flipping one spin amounts to a change 
in the local magnetic field experienced by another nearby spin) is replaced by an indirect 
interaction where the change in the magnetic field at the location of spin 2 is brought about 
by a change of its equilibrium position (caused by the spin flip of ion 1) in a spatially varying 
external field. 

For a more quantitative description we consider the force F — {{jl ■ V)i3) acting 
on an atomic magnetic moment jl — {'-^h/2)(j associated with spin {'h/2)(T in a magnetic 
field B. Without loss of generality we restrict the following discussion to the case where 
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Figure 2. Energy levels of two spins as a function of thei 
Si, i = 1,2 indicate the spin state. The resonance freque 
while when they interact these resonances are shifted by . 
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B — Bq + bz with constant offset field Bq and gradient b along the z— axis such that 
the force F, = ((7?i/2)6cr(^)) chan ges sign upon flipping the atomic spin (() indicates the 
expectation value). Using the position dependent spin resonance frequency oj = •jbz, the force 
can be expressed as a function of the change of uj when moving the ion along the z— axis: 

= {h/2)d,uj{a^^)) . (8) 

For an ion of mass m confined in a harmonic potential, characterized by angular 
frequency v, this linear force that acts in addition to the trapping potential and Coulomb 
force, shifts the equilibrium position of the ion, upon flipping its spin, by an amount 
dz = ^Fzl {mv'^) from its initial equilibrium position (compare figure[3j. Thus, in a classical 
picture, the ion finds itself on the slope of a harmonic oscillator potential, instead of at the 
potential minimum, after its spin was flipped. Consequently, it will start to oscillate around 
the new equilibrium position. In this way the ion's internal dynamics and its motional degrees 
of freedom can be coupled, even if the radiation causing the spin flip doesn't impart enough 
linear momentum to excite the ion's motion. This coupling, quantum mechanically described 
by an effective Lamb-Dicke parameter rye// = dz//Sz [Mintert& Wunderlich2001,Mintert& 
Wunderlich 2003, Wunderlich 2002, Wunderlich & Balzer 2003], has recently been observed 
experimentally [Johanning et al. 2009]. 




Figure 3. A single ion confined in a potential that is the sum of a harmonic oscillator and a 
linear term. The sign of the linear term depends on the internal spin state (± ^ -Ei ) of the ion. 
Upon flipping the ion's spin its motion can be excited. 

In order to see how the desired J— coupling arises, we consider how an ion's shift of its 
equilibrium position affects its internal energy. This change of its internal energy is given by 

hJ = -Fzdz = -F2/(mj/2) cx {b/vf. (9) 

As a consequence of one ion's shift of equilibrium position, through Coulomb interaction, 
the position and energy of neighbouring ions is shifted, too, giving us the desired (indirect) 
spin-spin interaction (compare figure|4|. 

I Such a field may be approximately realized, even though it does not fulfill V • B = 0, for example, on the axis of 
rotational symmefiy of a quadrupole field sufficiently far from the point of inflection symmetry. 




Figure 4. Sketch of two ions confined in a harmonic potential with superimposed spatially 
varying magnetic field B to illustrate spin-spin coupling mediated by Coulomb repulsion. 
When flipping the spin of ion 1, the state-dependent force acting on it will change sign, and 
thus, its equilibrium position will change (1 — * !')■ Through Coulomb interaction the second 
ion, too, will change its equilibrium position (2 — > 2'), and, since it moves in a field gradient, 
consequently will change its internal resonance frequency from uJ2 to ojg ■ 



Calculation of spin-spin coupling After having introduced this general physical mechanism 
that may be used to induce an effective spin-spin interaction, we will now briefly outline how 
this may be employed to engineer a variety of interactions between a collection of trapped 
ions. We consider the internal energies of N ions all of which experience a shift of their 
internal resonance frequency that, for now, depends on their axial position denoted by z, only: 

^ N N 

= - 2 ^ '^"'"^"^ " 2 ^ '^"'"^"^ ^^^^ 

n—1 n— 1 

Here, (?„ is the deviation of ion n from its equilibrium position Zn^ , 



d^UJn = dzUJ 



(11) 



and the expansion of u;„ has been included up to first orderly. 

The external (motional) degrees of freedom of N ions are determined by the applied 
trapping potential and their mutual Coulomb repulsion. Expanding the total potential up to 
second order gives 

3 N 

y(harm) = ^ ^ ^ A^a,.,j Qa. Qa.j (12) 
a — 1 — 1 

where the three spatial degrees of freedom a = x,y,z are uncoupled [James 1998]. 
The Hessian matrices Aaa,ij that characterize the potential experienced by the ions are 
real valued and symmetric in i and j and are transformed into diagonal matrices Da via 
Da = Sa^ Aa^ Sa with eigenvectors Q. 

In terms of the normal coordinates we obtain the following Hamiltonian describing 3A^ 



§ We assume a constant field gradient, that is, second and higher order derivatives of the field vanish. However, 
the expression [To]is still vaUd even if these derivatives are different from zero as long as the gradient doesn't vary 
appreciably over the extent Az of an ions spatial wavefunction. 
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uncoupled harmonic oscillators [James 1998] 

3 



a=l 



N N 



(13) 



where lya.n is the angular frequency of the normal vibrational mode n in direction a ~ x,y, z 
and the relation between the normal coordinates Qa,n and the local coordinates „ is Qa = 
{Qa,i, Qa.N) = Qa with = {qa.i, fc.w), furthermore Pa.n = m Q^.n- 

The total Hamiltonian includes internal and motional parts H = Hj + H^^\ where, for 
now, we restrict our considerations to the spatial z— direction where a field gradient is present. 
If a magnetic field gradient was present also in the radial direction (which is the case, e.g., for 
a quadrupole field close to its point of symmetry), then the radial vibrational modes would 
also contribute to the J— coupling. However, this contribution usually is small, since in order 
to maintain a linear ion string, the angular frequencies ly^.i and Uy i, characterizing the radial 
trapping potential, are usually chosen much larger than i^z.!- Taking into account that the 
J— coupling is proportional to the inverse square of the trap frequency (taking equation|9]as a 
guideline), we neglect the contribution from transverse vibrational modes for now. In [Porras 
& Cirac 2004a, Deng et al. 2005] transverse modes are explicitly considered (see below). 

After expressing the local coordinates qn in equation [TO] by the normal coordinates, 
9" = Sill SniQi, noting that Qi = Azi{al + ai), Pi = imvi^A{a\ - ai), and performing 
a suitable unitary transformation of H, H = WHU, with 

[/ = exp I - J ^ Az,e„i(7(")p, | . (14) 



.1=1 



we obtain [Wunderlich 2002] 



n— 1 n—1 i>j=i 

where the vibrational motion and the internal states are now decoupled (ajj, a„ are the creation 
and annihilation operator of vibrational mode n). The spin-spin coupling constants are given 
by 

N 

Jij — ^ ^ ^n^in^jn (15) 
n=l 

where the dimensionless constants 

AZndzLUi 

S'in = i>in (lo) 

have been introduced. The numerator of the fraction on the right-hand-side of equation [16] 
gives the change of the internal energy (or frequency) that ion number i undergoes when its 
position is shifted by an amount equal to the extension of the ground state wavefunction, 

Azn — \J 2mi, of vibrational mode n. This is set in relation to the energy (frequency) of 
a phonon of mode n. The matrix element Sin gives the scaled deviation of ion i from its 
equilibrium position when vibrational mode n is excited. Thus, tells us how strongly ion 
number i couples to vibrational mode n when its spin is being flipped. The sum in equationfTsl 
extends over all vibrational modes, since the local deviation, q^ of ion i from its equilibrium 
position potentially (if Sin is different from zero) excites all vibrational modes of the ion 
string. 
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(17) 



Noting that Vn oc vi and F^^i oc d^iOi and plugging the expression for Az into equation 
[161 we find from equationfTSjthat 

F F ■ 
b_ 

in agreement with the result |9] that was found from a simple model based on classical 
considerations (here, i^i is the angular frequency of the centre-of-mass mode of N ions). 
Going from the first to the second line in equation[T7] we have assumed a force caused by a 
constant magnetic field gradient causing a linear Zeeman shift of the spin states of ions i and 
j- 

This spin-spin coupling is only weakly sensitive to thermal excitation of the ion string. 
Thus, usually simple Doppler cooling will be sufficient to avoid unwanted thermal effects on 
the coupling constants Jij [Loewen & Wunderlich 2004]. Detailed calculations on this aspect 
will be published elsewhere. 

In [Ospelkaus et al. 2008] it is shown how an oscillating magnetic field generated in a 
micro-structured surface trap may be used to induce spin-coupling between ions and how this 
could be employed to carry out quantum gates. 



Spin-spin interaction created by an optical force The J— coupling described above may 
also be created by a state-dependent optical force (instead of the magnetic force assumed 
above) induced by an off-resonant standing laser wave [Porras & Cirac 2004a, Deng 
et al. 2005]. Applying state-dependent forces (e.g., optical or magnetic forces) along different 
spatial directions, combined with various trapping conditions, a variety of different coupling 
Hamiltonians of the general form 

Hi'°"^ = -\ E Y.J^f''>?-J2''^^f (18) 

a—x,y,z iyj % 

may be realized [Porras & Cirac 2004fl, Deng et al. 2005]. 

In particular, in [Deng et al. 2005] it is shown how state-dependent forces that act in 
the radial direction of a 1-d ion string may be exploited to simulate a transverse quantum 
Ising model and an XY-model, either with Coulomb crystals confined by a global potential, 
or linear arrays of microtraps where each ion is confined in its own, externally controlled 
potential well. It is demonstrated (by means of numerical and analytical calculations) under 
which conditions quantum phase transitions should be observable with 1-d strings of trapped 
ions, and what the expected characteristic features of these transitions are, as compared to the 
usual Ising and XY-models where nearest-neighbour interaction only is considered. 



2.1.3. Engineering the spin-spin coupling A straight forward way to apply a magnetic 
gradient that induces the long-range J— coupling described above in a linear Coulomb 
crystal is the use of current carrying coils in an anti-Helmholtz type arrangement, mounted 
such that their axis of rotational symmetry coincides with the trap z— axis. The field 
generated by these coils is easily calculated analytically and gradients up to several tens 
T/m are attainable with coils with a diameter of typically 1 mm. Some examples for 
possible configurations and computed examples for J— coupling constants are presented 
in [Wunderlich & Balzer 2003, Wunderlich et al. 2005, Wunderlich & Wunderlich 2009]. 
Such coils may also be integrated into micro-structured ion traps that are currently being 
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developed. Here, larger gradients are realistic (several hundreds of T/m) [Briiser n.d.]. When 
using current-carrying structures for generating magnetic fields, the currents, and thus the 
J— coupling, can be turned on and off as required, however, the stability of the magnetic field 
needs attention to avoid dephasing of spins. Alternatively, permanent magnets may be used 
[Johaiming et al. 2009]. 

Global trap potential If the ions are confined in one common trapping potential, then the 
strength of the J— coupling may be controlled globally, that is, all coupling constants may be 
varied together by adjusting the gradient b and/or the secular trap frequency Vi. In a linear 
Paul trap the trap frequency along the z— axis of a linear trap is adjustable by applying a 
suitable voltage to the end cap electrodes, while the radial trapping frequency can be adjusted 
by varying the amplitude of the applied rf trapping field. If the gradient h is generated by 
electrical currents, then these currents can be varied, and turning the magnetic field on and 
off using a desired envelope allows for simultaneous switching of all coupling constants. If 
permanent magnets are used, and their position is fixed along the 2— axis (typically generating 
a quadrupole field), the ions may be moved collectively along the 2:— axis to change the 
gradient they experience. Moving the ions in a linear Paul trap is achieved by applying suitable 
voltages to the electrodes that are responsible for the axial confinement. Even an always- 
on interaction (that arises, for instance, when using permanent magnets to generate a field 
gradient) used in conjunction with coherent manipulation of individual spins maybe employed 
such that prescribed entangled states can be generated in a time-optimal way [Fisher n.d.]. 
for the axial confinement. Even an always-on interaction (that arises, for instance, when 
using permanent magnets to generate a field gradient) used in conjunction with coherent 
manipulation of individual spins maybe employed such that prescribed entangled states can 
be generated in a time-optimal way [Fisher n.d.]. 

Locally variable trap potential The advent of micro-structured traps brings with it the 
possibility to attain control over local trapping potentials and thus, of engineering the 
J— coupling locally. For example, employing a linear Paul trap with dc electrodes that are 
segmented such that individual electrodes have an axial extension of the order of the inter-ion 
separation allows for shaping the local electrostatic potential such that long-range coupling is 
strongly suppressed and essentially only uniform nearest-neighbour coupling exists [Hugh & 
Twamley2005]. 

Furthermore, micro- structured traps with segmented electrodes allow for detailed local 
sculpting of J-coupling constants. This can be useful to generate effective time evolution 
operators making it possible to reach a desired entangled state in a single time evolution step 
without the need for refocussing. Examples how to generate cluster states in this way are 
given in [WunderUch & Wunderlich 2009]. 

A possible drawback of the use of small electrodes is that this usually implies an equally 
small distance between ions and nearby surfaces (i.e., typically of order 10 pm) which leads 
to high heating rates of the ions' secular motion, and thus limits the useful time available 
for coherent dynamics [Turchette et al. 2000]. Using ion traps embedded in a cryogenic 
environment considerably reduces this heating rate [Deslauriers et al. 2006, Labaziewicz 
et al. 2008], and therefore, may provide a feasible approach to using small electrode structures. 
A micro-fabricated linear surface trap cooled to 4 K that includes current-carrying structures 
for generating a magnetic field gradient was recently demonstrated [Wang et al. 2008]. The 
magnetic gradient was used to address individual ^*Sr+ ions in frequency space on an optical 
transition between Zeeman states of the electronic ground state 5Si/2 and the metastable state 
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4D5/2, respectively. 

Segmented ion traps also allow for transporting ions between different trap locations 
[Wineland, Monroe, Itano, Leibfried, King & Meekhof 1998, Kielpinski et al. 2002, Rowe 
et al. 2002, Home & Steane 2006, Hensinger et al. 2006, Schulz et al. 2006, Reichle 
et al. 2006, Hucul et al. 2008, Huber et al. 2008]. For moving ions, the required electrode 
structures may be considerably larger, typically of order 100 /im, thus helping to reduce 
unwanted heating. The transport of ions over small distances in segmented ion traps gives 
a further handle to engineer effective couplings by letting only a subset of ions interact at a 
time in a prescribed way. This could be used, for example, to generate two-dimensional cluster 
states using a one-dimensional linear Coulomb crystal [Wunderlich & Wunderlich 2009]. 

A possible architecture for a two-dimensional surface trap array exploiting magnetic 
field-gradient induced spin-spin coupling as described above was proposed in [Chiaverini & 
Lybarger 2008]. In [Clark et al. 2008] confinement of ^®Sr+ ions in a 2-d ion trap array with 
individual sites spaced apart by about 1.6 mm is reported. Even though a trap array with 
such a large inter-ion spacing does not allow for a useful spin coupling between the ions, 
it is nevertheless useful as a proof-of-principle and makes possible the investigation of the 
scaling behaviour of this trap array. It turns out that scaling down of the layered rf lattice 
trap alone is unlikely to be a promising route to achieve useful J— coupling. It appears that 
further modifications of this trap design are useful to ensure low enough secular frequencies 
for sufficiently strong J— coupling (compare equationfTTIl. 

Recently, two-dimensional and linear arrays of Penning traps that confine individual 
electrons have been proposed as a promising physical system for quantum computation and 
quantum simulations. Here, the physical origin of the J— coupling again is a force that 
depends on an electron's spin state and is caused by a magnetic field gradient and, as in 
ion traps, the Coulomb force between electrons mediates the spin-spin coupling [Ciaramicoli 
et al. 2005, Ciaramicoli et al. 2007, Ciaramicoli et al. 2008]. Thus, various couplings as in 
equation[3]with variable range could be realized with trapped electrons and used, for instance, 
to investigate quantum phase transitions. 

2.7.4. Simulating an effective transverse field An effective transverse magnetic field can be 
simulated by driving transitions between the spin states of ion i employing a radio frequency 
or optical field resonant with this ion's qubit transition [jj]; 

Ht= - hQRCOsiJ'h + (/)')cri'^ (19) 

^ _ (a^ + a^^) (e'"'"'*+^' + 6"^'^'"'+*') . (20) 

Here, flu is the Rabi frequency that for a linearly polarized rf field cos{uJxt + (f>') reads 
iljj = and the resonance condition is oj^ = w^*^ 

First transforming Ht by U, as was done above for the Hamiltonian H (equationflSll, then 
carrying out a transformation into an interaction picture with respect to — | J2n=i ^nO'z^'' + 
% i^noliO'n [Wunderlich 2002], expanding the resulting Hamiltonian up to first order in 

and a, neglecting terms that contain Sin 0, and applying the rotating wave approximation, 
one obtains for spin i 

iJ« = - ^^R{cos(^a^^ - sin0a«) (21) 

II Optical fields may be used to drive transitions in the optical range, but also in the RF range by applying a two- 
photon Raman excitation 

% In [Deng et al. 2005] the effects of residual spin-phonon coupling caused by the effective transverse field ai'e 
calculated 
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when the rf field or optical field is on resonance. Allowing for a detuning oj^f — i^^^^ between 
the applied field and the spin flip resonance frequency introduces a term proportional to 
(wr/ — w^'')cr'^'' on therhs of equationlSTI Thus, by driving spin transitions of individual ions, 

a variety of additional local effective magnetic fields B^i cx (cos(/!)„, — sin(/)„, [uj^f — oj^'^^Y' 
in Hamiltonian|3]may be simulated. In particular, by choosing the phase 0' of the rf field such 
that = and fulfilling the resonance condition for all ions, a global transverse field B'^ for 
the transverse Ising model|4]is implemented. 

2.7.5. Spin-boson model The spin-boson model ([Leggett et al. 1987, Weiss 1999] and 
references therein) is the starting point for many theoretical investigations modeling, for 
instance, the decoherence of a variety of physical system. It describes the dynamics of a spin- 
1/2 coupled to a bath of harmonic oscillators that represent the spin's environment. Thus, 
the spin-boson Hamiltonian contains a term proportional to CTz (describing the spin's energy), 
a term proportional to ax (describing spin flips), a term describing the harmonic oscillators, 
?i vi {a\ai + 1/2), and finally a term coupling the spin to the oscillators 

|a,^QQ; (22) 

where Qi is the position coordinate of the n— th harmonic oscillator while C/ describes the 
coupling strength between this oscillator and the spin. We note that all these ingredients are 
already present in the formulae presented above: Picking out a particular ion i, equation [TOl 
may be rewritten as 

H\=- - E a.-. g„ 

n— 1 

= f -.-i^^ - ^-i'^ E 9^^^ E (23) 

1=1 \n=l I 

and we can set C; = 9zWi {^2!^=\ Sni^ to obtain in equation l22l the coupling between spin 
i and the harmonic oscillators represented by the collective vibrational modes of a string of 
N ions. The term proportional to is obtained by driving transitions between the two spin 
orientations (see below, eq. 1211 1. 

Now the question arises how to physically single out a particular spin i from the sum 
over all spins in the second term of Hj (equationfTOb such that we are left with H} in equation 
I23] This could be done by exerting a state-dependent force only on ion i as is shown in [Porras 
et al. 2008]. Porras et al. detail how laser beam(s) focussed onto one particular ion give rise to 
optical forces that induce the required coupling. Furthermore, it is shown how phonon baths 
characterized by a variety of spectral densities may be simulated and how effects caused by 
the finite size of the ion crystal may be observed. 

2.2. Phase transition 

2.2.1. Quantum transverse Ising model The quantum transverse Ising model (TIM, equation 
m displays a phase transition already in one dimension [Sachdev 1999], even at zero 
temperature, between an ordered state and a disordered state upon varying B'^/ J. Milburn 
et al. [Milburn 1999, Milburn et al. 2000] proposed an implementation of a unitary map 
closely related to the TIM using Raman laser pulses. In [Porras & Cirac 2004fl] it is proposed 
how to observe this quantum phase transition (QPT) with trapped ions interacting with a 
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standing laser wave that creates a long-ranged spin-coupling as outlined above (section |2.1.2| i. 
Barjaktarevic et al. [Barjaktarevic et al. 2005] proposed to implement a unitary map that 
displays a QPT belonging to the same universality class as the transverse Ising model. 

In order to get a qualitative picture of such a phase transition we first consider the 
collection of spins described by Hti (equation IU in the limit where the effective transverse 
field, assumes a finite value, but the mutual coupling ,/ between spins vanishes, B'^ ^ 
= J, such that, at zero temperature, all spins are oriented along the x— direction and the 
groundstate |0) — |— >)j — ni(IT)i + of th^ system is a product of eigenstates of 
erf. Consequently, (0| cricr^j |0) — Sij, that is, the orientation of different spins with respect 
to the the z— axis is uncorrected. On the other hand, if the spin-spin coupling doesn't vanish, 
B'^ ^ \J\ > 0, the spins will acquire a tendency to align along the z— axis. Thus, the 
ground state of the system will start to exhibit non-vanishing, but short ranged correlations 
between spins located at sites and rj, respectively, with a characteristic length scale £, 
[Sachdev 1999] 



Here, |0) is the system's ground state when ^ | J| ^ 0. 

Now, we turn to the limit B'^ ^ J > 0. If B'^ vanishes completely, then an infinite 
system of spins possesses two degenerate ferromagnetic ground states where all spins are 
either parallel or anti-parallel to the z— axis, and the total magnetization of the system is Mq. 
Turning on a small field B'^ <^ J pointing in the x— direction will tend to reorient some of the 
spins along the x— direction (classically, their magnetic moments will start to precess around 
the direction, thus flipping spins between the two possible z— orientations). The ground 
state of the system, |0) is still close to the B'^ — case and is obtained from a perturbative 
expansion in the small parameter B'^/J. One finds (0| erf |0) = M/Mq where M < Mq is 
now the spontaneous magnetization of the system. Again calculating the correlation of the 
orientation of different spins in the z— direction for this ground state gives [Sachdev 1999] 



The expressions l24l and [25] for the spin correlation are incompatible in the sense that a 
continuous variation of B'^J J between the two limiting cases considered above doesn't allow 
for an analytical connection between the two expressions on the right-hand-sides of equations 
l24l and |25] Thus, there must exist a critical value of B'^j J where the correlation changes 
between these two expressions. This is the case even at zero temperature being indicative for 
a quantum phase transition. 

An experimental realization of the proposal by Porras and Cirac [Porras & Cirac 2004a] 
with two trapped and laser cooled ^^Mg"*" ions was recently reported by Friedenauer et al. 
[Friedenauer et al. 2008]. Here, the coupling J between spins is induced by a state-dependent 
ac Stark shift caused by a "walking" standing wave. Two hyperfine states of the electronic 
ground state of each ^^Mg+ ion served as an effective spin-1/2 with eigenstates ||) and \\). 
In brief, the experimental procedure was as follows: Optical pumping is employed to first 
initialize the two ions in state Resonantly driving the transition between the two 

hyperfine states with RF radiation simulates an effective transverse field Blj, (section 12.1.41 ) 
and allows for preparation of the initial paramagnetic state ^) (both spins are oriented 
along the effective B'^ field). Then, the J— coupling between the two spins is ramped up 
adiabatically from zero to a desired value such that the ratio J jBlj. < 5.2. For J ^ B'^ one 
expects that the two spins are found in a ferromagnetic state, that is, aligned along the z— axis. 
In general, the new ground state will be a superposition of the states \i, i) and |t, T) which 



(0| cTf ct| |0) cx e 



ri-rj\/i 



(24) 




(25) 
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is indeed found experimentally. In a further experiment the two spins were prepared initially 
anti-parallel to the effective B'^ field, an excited paramagnetic state. Turning on the spin-spin 
coupling adiabatically then resulted in an excited state in the presence of strong ./—coupling 
{J jB'^ » 1): the anti-ferromagnetic state |t,i) + |i,T)- In both cases the entanglement of 
the final state was proven experimentally. 

2.2.2. Scalar fields Classical phase transitions between spatial configurations of ions 
confined in a linear or ring shaped trap have been investigated experimentally and theoretically 
[Waki et al. 1992, Raizen et al. 1992, Birkl et al. 1992, Dubin 1993, Schiffer 1993, Drewsen 
et al. 1998, Kjsergaard & Drewsen 2003, Morigi & Fishman 2004fl, Morigi & Fishman 2004b, 
Fishman et al. 2008]. Retzker et al. [Retzker et al. 2008] suggest to explore a phase transition 
between a linear and a zig-zag configuration [Fishman et al. 2008] in the quantum regime, that 
is, with the ion string cooled close to its vibrational ground state. The ion trap parameters are 
chosen such that a linear Coulomb crystal forms (i.e., the confinement in the axial direction is 
weaker than in the radial direction), and, in addition, i^y i 3> I'x.i such that one radial degree 
of freedom is essentially frozen out after laser cooling, but the other may be employed for 
quantum simulations. Linear Paul traps as well as Penning traps where the ions are confined 
on the axis of rotational symmetry [Powell et al. 2002] are well suited for this purpose. 

For //l a; below a critical value i^f ^, a zig-zag configuration of the ions will form with 
two degenerate states. It is shown in [Retzker et al. 2008] that coherent tunneling between 
these two states maybe observed and used for interferometry. Furthermore, by adjusting lyi^x 
around ^ the parameters of a Hamiltonian describing the ions' motional excitation along 
the a;— direction may be determined. This Hamiltonian describes at the same time a non- 
linear scalar field and by adjusting i^i ^ the effective mass of this field may be varied, in 
particular in the vicinity of a phase transition. An attractive feature of this proposal is that 
global adjustment of the radial trapping frequency i/i^x is sufficient and no local control over 
the trapping potential is needed. 

2.3. Neural Network 

In [Pons et al. 2007, Braungardt et al. 2007] it was shown that a linear Coulomb crystal of 
laser cooled trapped ions may be used to implement a neural network (NN) - according to 
the Hopfield model [Hopfield 1982] - allowing for robust storage of classical information, 
and furthermore, that such a NN may be employed for error-resistant quantum computation. 
We will first briefly introduce some features of the Hopfield model (see, e.g., [Rojas 1996] 
chapter 13 and references therein) and then outline how some features may be realized with 
trapped ions. 

A neural network according to the Hopfield model consists of N nodes (corresponding 
to neurons) that are all mutually interconnected and we indicate the (synaptic) strength of the 
connection between node i and node j by . This connection is symmetric, that is = J^j 
and the nodes do not feed back to themselves, Ja = 0. Below we will see what exactly the 
meaning of the coupling strengths is. Each node can assume one of two possible states 
which we label by = ±1 (compare figure|5]l. When modeling neurons, this would indicate 
whether the neuron is "firing" (a sudden change of a biologically relevant potential difference) 
or in a quiescent state at a given time. Furthermore each node is characterized by a threshold 
hi that, in addition to the couplings Jy , influences its response to the the input from other 
nodes. 

A characteristic feature of a Hopfield NN is that it can store information distributed 
over all nodes. Thus, it is robust against errors in individual nodes or even the loss of 
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some nodes and against thermal noise. The dynamics of the NN, being determined by the 
strengths of interactions between all nodes, brings the system back to one out of a selection of 
particular states (that we call patterns) once this state was disturbed. This feature may be taken 
advantage of for pattern recognition: When the NN is presented with some input state, which 
is done by initializing all nodes with a given value, then the NN evolves into the pattern that is 
was trained to recognize that is closest to this initial state. So how does the NN dynamically 
evolve? 




Figure 5. Illustration of a neural network that consists of four nodes according to Hopfield 
[Hopfield 1982]. Each node inay exist in one out of two states indicated by Si = ±1, 
i = 1,2,3,4. The nodes are all interconnected with Jij = Jji. Here, the network forms 
a pattern g = (—1, +1, +1, —1)^- (Note that the number p of patterns that can be stored by 
a Hopfield network is about p = 0.14A'^ where N is the number of nodes). 



We start at time t with an initial global state of the NN labeled by the iV— component 
vector s where each local node is characterized by s;. Then the state of each node at time 
t + dtis determined by the sign of the activation function A at time t + St: 

+ =sgn(A(t + (5t)) , (26) 

where 

A{t + &)=Y^ AjSj {t) + . (27) 

This means that at time t + St, node i will be firing (si — +1) or in a quiescent state (si = — 1) 
depending on the state of all other nodes at time t, weighted with the strength J^ , and the 
individual threshold hi. One of the two following alternative model assumptions is usually 
made: The updating of node i may either occur sequentially, one node after the other, or in 
parallel all at the same time. 

We now define the cost function 

Hh{s) = -'^ ^ JijSiSj -"^hiSi (28) 

and it can be shown [Rojas 1996] that an evolution of the network according to the dynamical 
rule for sequential updating |26] (together with the definition |27l ) leads to a monotonous 
decrease of Hh This means that flipping spins (in agreement with the rules) drives the 
system towards local or global minima in the "landscape" defined by Hh over the space of s. 

In this brief account we consider the system at zero temperature only 
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Immediately we realize that the expression defining Hjj is formally equivalent to the spin 
Hamiltonian Hg (equation[3]) considering the case of spin- 1/2 particles with j\f = = j\f . 
As noted at the beginning of the section on spin models, further simplifications of Hs, namely 
uniform J— coupling, only nearest-neighbour interaction, and no transverse field lead to the 
Ising model (equation [ij. When considering the Ising model one usually assumes that all 
coupling constants are positive or all are negative and thus arrives at a limited set of minima 
(e.g., for J > all spins are all either up or down). However, here the couplings Jy are not 
restricted to being all positive or all negative, thus allowing for a rich "energy landscape" with 
many minima. 

Our goal is to "teach" the NN a specific pattern, ^, i^j = ±1 towards which the NN 
is supposed to relax, if it is displaced from this pattern, or if presented with a similar input 
pattern. In other words, this pattern must correspond to a minimum of Hh- Where exactly 
this minimum occurs, that is, for which particular configuration s, is determined by the choice 
of the interaction strengths . Therefore, if one wants the NN to "remember" a particular 
pattern ^, (i.e., create a minimum of Hh for a specific arrangement ^ of the nodes), the 
J— couplings need to be adjusted accordingly. At first sight (and probably still at second) it 
is not obvious how a set of J— couplings may be found that creates a minimum of Hh for 
a desired pattern and thus makes the NN evolve towards this pattern. It turns out that using 
Hebb's rule [Hebb 1949, Mezard et al. 1987] to determine the values of the couplings 
gives just the right set of J— couplings. 

In order to memorize a pattern ^, each is modified according to Hebb's rule of 
learning, such that 

= A J°/^ + e^;^, . (29) 

When starting with an "empty" memory, then J"^''- — 0. In the original Hopfield model 
A = 1 = e; another choice is e = 1 /N . With the set of interconnections, determined by 
equation |29] the NN, upon turning on its dynamics, will eventually end up in the memorized 
pattern ^. If the NN is supposed to memorize p patterns, ^'^ with jx— 1,2,3, instead of 
just a single pattern, then the learning rule is modified, namely 

-/^.—^Ecrer (30) 

Here, we start with an empty memory, i.e. initially Jy = OVi, j. After fixing the strength 
of the pairwise J— coupling in this way, the patterns are "memorized" by the NN, and, when 
presented with some initial state, the NN will evolve towards a minimum of Hh associated 
with one of these patterns. The upper limit for the number of patterns that such a network 
may faithfully recall via its associative memory is about Q.IAN [Rojas 1996]. 

Now, it will be discussed how a linear Coulomb crystal of trapped ions may be employed 
to realize a Hopfield-type NN. We have already noted the similarity between Hh (equation 
|28] | on one hand and the two spin Hamiltonians i/i"' (equation |3]l and (equation [TsTi 

for jI^'^ = = Jj^J^ on the other. A further similarity can be seen between equation l30l 
defining the synaptic strengths in the Hopfield model and equation [15] where the ion spin 
coupling constants are given. We rewrite equation [TS] using the expression [T6l assuming a 
constant field gradient and obtain 

ji^n) ^n^^^y^^^ (31) 
n—l 

where the matrix elements Sm may assume positive and negative values, however usually 
different from ±1 as in the generic Hopfield model. Also, the sum in equation [3T| extends 
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over all N vibrational modes instead of p patterns, and the contribution to J^- ' from each 
vibrational mode n is weighted by the square of the respective vibrational frequency 
Despite these differences, it is possible to store spin patterns with trapped ions in a robust 
way. 

What are the patterns, based on the learning rule [31] that can be stored by a trapped ion 
neural network? The orientation of a specific ionic spin i in a pattern ^ of spin orientations is 
determined by the sign of Sin for a given vibrational mode n. In case of the COM mode, for 
example, we have Sin > (or Sin < 0) Vi which means that in the two associated patterns 
all spins are either pointing up or all down. 

Numerical simulations [Pons et al. 2007, Braungardt et al. 2007] show that, in order to 
robustly store a pattern of spin orientations it is necessary to adjust the trapping potential 
such that the vibrational frequencies vi and 1^2 (that appear in the learning rulelSTTi are nearly 
degenerate. This can be achieved by using ion traps with segmented electrodes where the 
Hessian of the trapping potential can be tuned by applying suitable voltages to individual 
segments [Hugh & Twamley 2005, Wunderlich & Wunderlich 2009]. The robustness of the 
ion NN has been quantified by Monte Carlo simulations for a string of 40 "^^Ca^ ions trapped 
such that the two lowest vibrational modes characterized by h'l and 1^2 are nearly degenerate. 
For this purpose, a spin pattern fj, = 1,2 corresponding to either one of the two modes is 
prepared as the initial condition of the ions' spin orientations. Then a fixed number r of spins 
are randomly selected and their state is inverted. This is followed by the time evolution of the 
ion string until an equilibrium state is reached. The probability of the ion string for ending up 
in the initial pattern is found to be more than 97 % even when up to eight randomly selected 
spins are initially in the wrong state (i.e., they are not matching the desired pattern). Thus, 
even with 20% of the spins being flipped, the ion string returns towards a desired pattern. 

So far we have looked at the robustness of a Hopfield-type neural network where each 
node consists of a single ion and which is used for classical information storage. For a 
quantum neural network, the elementary unit for information processing is, instead of a single 
spin- 1/2 with two eigenstates corresponding to |0) and a pattern xi of N spins. In [Pons 
et al. 2007, Braungardt et al. 2007] it is detailed how an appropriately chosen set of patterns 
with eight trapped ions may be used as qubits thus forming a quantum neural network. In 
particular, single-qubit gates and two-qubit gates by adiabatic passage are worked out that 
make use of additional effective fields. 

2.4. Mesoscopic systems 

The question if and to what extent quantum mechanical effects play a role in mesoscopic or 
even macroscopic systems was posed already shortly after quantum theory was developed 
[Schrodinger 1935]. Except concerning some particular physical phenomena, for example, 
macroscopic currents in a superconducting material, it is usually assumed that superposition 
states, interference, and entanglement are not relevant when dealing with entities made 
up from a large number of individual constituents that, in addition, may interact with a 
noisy environment at non-zero temperature, for example, biologically relevant molecules. 
Indications that quantum mechanics indeed plays an important role even for large molecules 
come from recent experiments where interference with C70 molecules was observed 
[Hackermuller et al. 2004]. 

A theoretical study of how dephasing caused by coupling to a noisy environment can 
actually assist in transferring energy (or information), instead of impeding it, across a network 
of quantum nodes is presented in [Plenio & Huelga 2008]. It is shown that, in a simplified 
model, light-harvesting molecules may take advantage of this process. An experimental 
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verification directly with these molecules is difficult at present. However, trapped ions may be 
used to implement a suitable model to test such a prediction. For a first simple demonstration, 
appropriately chosen internal levels of a single ion may be used to model the nodes of 
a network and induced and spontaneous transitions between these levels to mimic energy 
transport. 

Another recent theoretical study demonstrates that entanglement may persist in a model 
two-spin molecule even when it is exposed to a decohering hot environment, a situation 
typically encountered for biological molecules [Cai et al. 2008]. The time dependent 
Hamiltonian modeling the two-spin molecule used for this study is 

HM{i) = J{t)alal + B{t){al + al) . (32) 

This molecule is then coupled to a noisy environment modeled by a master equation. After 
performing a unitary transformation such that ax a^ and a^ —ax we recognize a time- 
dependent variant of the quantum transverse Ising model |4] which could be realized with an 
ion spin molecule by using the spin-coupling present in equation [15] and by adding a driving 
field as in equation 1211 

3. Quantum Optics Simulating Quantum Optics 

The formal similarities for creation and annihilation operators for photons and phonons on 
one side, and spin S* = 1/2 particles on the other side suggests the equivalence of various 
and seemingly different experiments within the field of quantum optics. In this notion, 
a Ramsey measurement might be regarded as simulating a Mach-Zehnder interferometer, 
since a wavefunction is coherently split and recombined after a independent time-evolution 
of its constituents. Therefore one general proposition which holds also for other classes 
of experiments discussed in this paper is maybe of particular importance here: simulation 
goes both ways. If the Hamiltonian describing a setup in an ion trap, possibly including 
static and radiation fields, is equivalent to the Hamiltonian of another system, then the ion 
trap experiment can be thought as a simulation of the other system of interest. But on the 
other hand it also means, that these other system might reveal some insights about certain 
aspects and experiments in ion traps, and, depending on the specific technical difficulties 
when trying to perform it in an ion trap, it might also be imaginable to simulate some aspects 
of ion trap physics in other systems. As pointed out earlier, however, trapped laser-cooled 
ions are a flexible, well controlled and well understood system, which can be brought into 
a variety of regimes, where the physics resembles other systems. Nowadays we are familiar 
with the fact that both photons but also the internal states of particles can be entangled or 
combinations of them and it might sound unusual to phrase any of such experiments as a 
simulation of the other Nevertheless the creation and detection of entanglement of particles or 
phonons is subject to different problems and certain aspects of entanglement and its meaning 
in investigating fundamental questions of quantum mechanics might best be looked at in either 
approach [Wineland, Monroe, Itano, King, Leibfried, Myatt & Wood 1998]. 

3.1. Theoretical Background 

Already in the early days of quantum information it became clear that quantum computing at 
a level of factorizing large numbers would be a daunting task (see [Wineland, Monroe, Itano, 
King, Leibfried, Myatt & Wood 1998]), though, since then, astonishing progress was made 
this field. The first articles considering quantum simulations were, to some extent, motivated 
by the insight that the requirements for quantum simulations (in terms of gate fidelity etc.) are 
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less stringent despite their potential to investigate systems that are beyond reach of simulations 
with classical computers. Despite all progress, this statement still holds and is one motivation 
for this topical review. 

The experiments in quantum optics that are considered in [Wineland, Monroe, Itano, 
Leibfried, King & Meekhof 1998, Leibfried et al. 1997] can be simulated by trapped ions 
interacting with radiation fields that drive Raman transitions between hyperfine levels. A 
similar description can be obtained when this transitions are directly driven using microwaves 
(using auxiliary inhomogeneous magnetic fields) [Mintert & Wunderlich 2001, Mintert & 
Wunderhch 2003, Wunderlich 2002, Johanning et al. 2009] or when optical transitions 
between metastable states are used. For a harmonically bound atom interacting with two 
travelling wave light fields, both detuned from any real state and with frequency difference 
Aw and phase difference (f) the Hamiltonian in the interaction picture can be brought into a 
form H^i (after rotating wave approximation and adiabatic eUmination of the near resonant 
excited states and in the Lamb-Dicke limit 77^ <|(n + a^)) <C 1) [Wineland, Monroe, Itano, 
Leibfried, King & Meekhof 1998] with integer I indicating changes in the motional state, and 
e e [0, 1] indicating changes in the internal state. This Hamiltonians are frequently used in 
experiments with laser cooled ions; to name a few, iJio is called the carrier transition, H12 is 
the second blue sideband, ifoi is called coherent drive, Hq2 is the squeeze drive. Nesting 
and concatenation operators of the set Hoi,Ho2,Ho3,Hio,Hii,H-i'2,Hi3 is sufficient to 
efficiently generate a wide range of Hamiltonians, as, for example the Hamiltonian for a spin 
s = 1/2 particle with mass in an arbitrary potential [Leibfried et al. 2002]. 

These chaining of Hamiltonians could be used for a variety of experiments, as, 
just to highlight a few, the phonon maser [Wallentowitz et al. 1996] and three phonon 
down conversion [Wineland, Monroe, Itano, King, Leibfried, Myatt & Wood 1998] and 
has been used to create motional cat staes [Monroe, Meekhof, Leibfried, King, Itano & 
Wineland 1996], coherent states [Meekhof et al. 1996fl, Meekhof et al. 1996/?], and to create 
a Hamiltonian similar to two photon excitation in cavity QED [Leibfried et al. 1997]. An 
example discussed in the following subsection is an an experimental analogy to a nonUnear 
two beam (e. g. Mach-Zehnder) interferometer as proposed in [Wineland, Monroe, Itano, 
King, Leibfried, Myatt & Wood 1998]. Using ions for the implementations allows for high 
state preparation and detection efficiency and eliminates the need for data postselection. 

All the experiments above used only a single ion and innumerable possibilities come into 
reach using multiple ions. The first demonstration of the CNOT gate [Cirac & Zoller 1995] 
was a step in this direction using, however, a single ion and its motional degree of freedom 
[Monroe, Meekhof, King, Itano & Wineland 1995]. This concept can be adapted to entangle 
different atoms when it is carried out after the state of one ion is mapped onto the centre 
of mass (CM) mode and afterwards mapped back [Schmidt-Kaler et al. 2003, Leibfried 
etal. 2003]. 



3.2. Experiments 

As a demonstration of a quantum simulation with a single ion, the Hamiltonians 
Hii, H12, -ffi3 are used in [Leibfried et al. 2002] to simulate a set of two nonlinear beam 
spUtters forming a nonlinear Mach-Zehnder. A single ^i3e+ ion provides the internal state | a) 
with 

|0)„ = \F=l,mF = -l) 

(33) 

= \F = 2,mF = -2) 
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and the motional state |n)j labelling the nth vibrational state. The experimental sequence is 
depicted in|6]and is as follows: the ion is cooled to the ground state of motion |1)^ |0)^. A 
Raman 7r/2 pulse applied on the appropriate sideband is used to create the state proportional to 
I l)a' l'-')6' + |0)a' !"■)&' ■ This is equivalent to a beamsplitter: analogous to coherently splitting 
light into two separate directions, the population is coherently split into two states which, in 
this case, differ in both internal and motional quantum number. Subsequently the potential 
of the endcap electrodes is changed yielding trap frequency change /S.Vz.i for a time t which 
introduces a phase shift n/S.Vz,it = n(f) when the ion is in the nth vibrational level. Finally 
a second 7r/2 pulse on the same sideband is applied which maps the phase difference to 
the population in the internal states |0)^ and which is measured by Doppler cooling 
and fluorescence detection. The absolute acquired phases are irrelevant as is the total arm 
length of a Mach-Zehnder interferometer: longer arms only tend to make it more susceptible 
to perturbations, however, the signal at the outputs is given by the differential phase shift. 
The sequence is repeated to improve signal to noise ratio and to extract the sinusoidal time 
evolution in the different states. 

The experiment can be interpreted as nonlinear Mach-Zehnder interferometer. A 
nonlinear beamsplitter can annihilate a photon in mode a while creating n photons in mode 
h (or vice versa). The annihilation operator a is replaced in this experiment by the atomic 
raising operator cr^ between the states \F = 2,mp — —2) and \F — l,mp — —1). These 
two operators are not equivalent but they act the same as long as the {|0)^ , \ subspace 
is not left, which is fulfilled as long as the input state is |0)^. The optical mode with 
lowering operator b is replaced by the motion along a selected direction and phonon number 
states |n)j,. A ?ith order nonlinear beamsplitter is then given as 



When this operator brackets an operator that causes a phase shift in one arm by ncf), the system 
resembles a Mach-Zehnder interferometer and the probabihty of detecting a phonon in output 
a" becomes 



This experiment requires a sideband transition in the nth sideband, which can be very 
inefficient with increasing n as the sideband Rabi frequencies scale with 77" and the 
preparation becomes slow. An faster alternative route exploiting an auxiliary state and 
first order sidebands is detailed in [Wineland, Monroe, Itano, King, Leibfried, Myatt & 
Wood 1998]. 

An increase in fringe frequency proportional to the order of the nonlinearity is found 
but the fringe contrast is reduced with increasing n, as preparation and manipulation 
imperfections spoil the fidelity of the states and the detection efficiency is below one for all 
experiments. The sensitivity of the output signal to phase changes is maximized at the steepest 
slope of the signal and the signal to noise ratio for the n — 2 and n — 3 measurements is found 
to be below a perfect first order interferometer with unit detection efficiency. A potential 
improvement by increasing n from 2 to 3 is almost perfectly cancelled by the reduction of 
fringe contrast in the data presented in [Leibfried et al. 2002]. 

4. Transverse standing wave: Bose-Hubbard-physics (and beyond) 



n 



(34) 



Pa" - 



- [1 - COS {ncf)] 



(35) 



Laser cooling and trapping has been considered [Hansch & Schawlow 1975, Wineland & 
Dehmelt 1975] and successfully applied (see [Ghosh 1995, Metcalf & van der Straten 2001], 
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Figure 6. Ion trap Mach-Zehnder interferometer: the manipulation of the wavefunction of 
a trapped ion can be thought analogous to a set of two nonlinear beamsplitters forming a 
nonlinear Mach-Zehnder interferometer. 

and references therein) to both neutral atoms and ions. Both types of laser cooling experiments 
are often very much alike in terms of the description of the experiment (trap frequencies and 
depths, Rabi frequencies etc.) but also in terms of experimental techniques. A recent topic 
in cold neutral atom physics is the combination of cold trapped atoms and optical lattices 
which allows to experimentally investigate the Bose Hubbard model [Fisher et al. 1989], 
which is related to the Hubbard model in soUd state physics [Hubbard 1963]. The freedom 
in changing tunnelling and interaction parameters, the defect free optical lattices and the 
controlled introduction of disorder [Schulte et al. 2006] allow nice and clean experiments 
[Greiner et al. 2002, Morsch & Oberthaler 2006]. Also fermionic atoms have been loaded into 
traps and lattices and brought to quantum degeneracy. The description of such experiments 
by the Bardeen-Cooper-Schrieffer (BCS) theory [Bardeen et al. 1957, Giorgini et al. 2008] is 
known from superconductivity and this analogy encourages the hope that this experiments can 
help to investigate and understand high Tc superconductivity [Hofstetter et al. 2002, Holland 
et al. 2001]. 

The large and long range Coulomb force being the dominant force in ion traps has a 
number of consequences: trap frequencies, trap depths, interaction energies, and particle 
separations are typically higher by orders of magnitude compared to traps for neutral atoms. 
Thus Bose-Einstein condensation (and also Bose-Hubbard physics), one long sought holy 
grail in quantum optics, has therefore been investigated with neutral atoms but not with 
ions: when cooling neutral atoms at sufficient density, their wave functions can overlap and a 
macroscopic population of the trap ground state forms [Anderson et al. 1995, Davis, Mewes, 
Andrews, van Druten, Durfee, Kurn & Ketterle 1995]. When on the other hand an ion chain 
is cooled to the lowest vibrational mode [Diedrich et al. 1989], the ion separation is orders of 
magnitude larger than the spatial extension of the wave function. 

But quantum degeneracy and Bose Hubbard physics can still be simulated in ion traps 
when the ions take the role of lattice sites which can be populated by phonons, more 
specifically excitations of the vibrational modes transverse to the weak trap axis. The phonons 
therefore represent the neutral atoms located in one lattice site and phonons can tunnel 
between sites as do the atoms in optical lattices. Furthermore, as we will see later on, the 
phonons can be made interact by anharmonicities in trapping potential. This anharmonicities 
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can be introduced by the light shift of an off-resonant standing wave. This could be used 
to demonstrate a Bose-Einstein condensate (BEC) of phonons [Porras & Cirac 2004/?], the 
superfluid-Mott insulator transition [Deng et al. 2008] and, at higher levels of interaction 
energy, the Tonks-Girardeau gas [Deng et al. 2008] and even frustrated XF-models [Schmied 
et al. 2008]. In ID, BEC at non-zero temperature is only possible when the number of atoms 
(or here phonons) is finite, and we are in the weak coupUng regime, which means, that the 
average interaction energy is smaller than the expectation value of the kinetic energy. A 
peculiarity of the effective phonon interaction being dependent on the position of the standing 
wave is that it can be tuned, even in sign, similar as a the magnetic field can change the 
scattering length near a Feshbach resonance. 

Besides the general appeal to demonstrate Bose-Einstein condensation in yet another 
system, Bose-Hubbard physics in ion traps offers advantages, and also drawbacks, over 
present neutral atom experiments: experiments in ion traps do not only offer global control 
of tunnelling rate and interaction energy (as we will see later), equivalent to ramping the 
optical lattice depth or applying a magnetic field to obtain a certain scattering length near a 
Feshbach resonance. Ion traps allow also for local control: irregular lattices correspond to 
site dependent tunnelling rates and can be obtained in segmented microtraps; the interaction 
energy can as well be made site dependent by having independent standing wave intensities 
and phases at the different ion positions - in terms of neutral atoms this would correspond to 
arbitrary 3d Feshbach magnetic fields. Other interesting features include single site control 
and readout of the state and the potential to extract the density matrix either partially or, even 
fully (in small systems). On the other hand the experiments discussed here are restricted to 
bosonic phonons and thus Bose-Hubbard physics and today's ion traps are mostly restricted 
to one dimension but true two-dimensional systems can be obtained in surface traps [Seidelin 
etal. 2006]. 

4.1. Background 

When we speak of phonons in the context of Unear ion traps, we often think in terms of 
axial normal modes. In this situation momentum is transferred by local interaction with the 
addressed ion which couples efficiently to the motion of the whole chain due to changing 
separations and acting Coulomb forces. In this chapter radial modes of the ion motion are 
considered. Classically a radial displacement of a single ion excites both radial motion of 
other chain members but also the axial motion of the whole chain. But when the axial ion 
separations are much larger than radial excursion of the ion motion, the coupling of this 
motion to all other modes is very inefficient: the motion is almost perfectly perpendicular 
to the Coulomb interaction between ions all the time and does not change the ion separations 
to first order. Thus radial phonons can be viewed as being located to one ion but the non 
vanishing coupling makes the localization non permanent and re-emerges as a finite phonon 
tunnelling rate. 

The Hamiltonian describing the system reads [Deng et al. 2008] 



where the trapping potential Vt could be simply quadratic like the effective potential created 
by a simple linear Paul trap, or be array of microtraps achieved by micro structuring [Rowe 
et al. 2002, Stick et al. 2006]. We will look in more detail into one radial coordinate x 
and excitations at the trap frequency Vx,i in this direction. The trap frequency Vy^i can be 
chosen to be non-degenerate to allow for addressing specific sidebands and to avoid unwanted 




(36) 
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a) 




b) 





Figure 7. Experimental setup (not to scale) for the simulation of Bose-Hubbard physics. 
Images a) and b) are sketches of the trap and the applied standing wave potential. The 
illustrations are not drawn to scale, especially the wavelength of the standing wave is greatly 
exaggerated for visual clarity. Figure a) shows the trap as seen from a one of the endcap 
electrodes. Figure b) shows the setup from the side, as seen in most cases from an imaging 
detector; the weak trap axis is perpendicular to the standing wave. 



couplings. The absolute value of the Coulomb force between two particles separated by a large 
distance 5z in z-direction and displaced only slightly by 5x in the a;-direction is 



\F\ = 
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The component in .x-direction of the Coulomb force is given by 
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When the radial separation 5x is identified with Xi ~ Xj and the axial separation 5z with the 
unperturbed z° — 2:°, the radial potential for small relative displacements becomes 
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Introducing the quantities 



1 1 



(40) 
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the second quantized form becomes [Deng et al. 2008] 

The quantities tij are site dependent hopping terms resulting from the non vanishing 
coupling of the radial motion between all ions in the chain. Coupling of ion motions and 
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Figure 8. (colour online) The transverse potential as a sum of a parabolic trapping potential 
plus a far off-resonant standing wave: the black curve shows the unperturbed effective trap 
potential, while the red is the total effective potential with a standing wave having a maximum 
at the ion position, resulting in a effective repulsive interaction of phonons. The blue shows 
the total effective potential with a standing wave having a node at the ion position, resulting in 
a effective attractive interaction of phonons. 



therefore tunneling of transverse phonons is the stronger the closer the ions are. The next- 
neighbour tunnelling is therefore the strongest but its slowly decaying behaviour with lattice 
site difference can already be seen from equation ( [39l ). leading to long range tunnelling in 
stark contrast to pure next-neighbour tunnelling, as found in neutral atom experiments. The 
Coulomb repulsion from adjacent ions has a maximum for 5x — and counteracts the radial 
trapping potential, and the v^. \ are the reduced radial trapping frequencies. These changes 
in the radial trapping frequencies and therefore the reductions in the phonon energies are the 
strongest for small ion separations and thus lead to an axial phonon confinement towards 
the centre of the trap [Deng et al. 2008]. This confinement becomes the stronger, the more 
inhomogeneous the ion separations are. Therefore this effect is more pronounced in a linear 
ion trap compared to the isospaced array of microtraps (with almost constant ion separations) 
and becomes more evident with increasing ion number. 

Due to the large separation of ions and the small transverse oscillation amplitudes, the 
effect of the Coulomb interaction to the radial potential is just a change in the trapping 
frequency and the anharmonicities, that should in principle be already present, have been 
neglected so far. A purely harmonic trapping potential results in equally spaced oscillator 
levels, meaning that the energy of having two phonons located on different sites or both at the 
same site gives the same total energy. Interacting phonons means that the energy of having 
two phonons at the same site is the sum of the the two single phonon energies, altered by an 
interaction energy, and the oscillator levels are no longer equidistant. A nifty idea to mimic 
this interaction is to introduce an anharmonicity to the trapping potential and thus to break 
the regular spacing of the trap levels. This can in principle be done by any non harmonic 
potential as, for example a sinusoidal potential as known from optical lattices generated by an 



Quantum Simulations with Cold Trapped Ions 



26 



off-resonant standing wave with wave vector /cgvv [Deng et al. 2008]: 

JV 

^^sw = ^ I] cos2 [ksw X, + (42) 
1=1 

where F is the peak AC Stark energy shift. For small anharmonicities the Hamiltonian 
Hxo + Hsw can be expanded in the Lamb-Dicke-parameter rix,sw with 



ilxsw = kt.^/2mhvx,i (43) 

and is sufficiently approximated when truncated after second order terms and the Hamiltonian 
can be brought into the well known form of the Bose-Hubbard model [Deng et al. 2008, Fisher 
etal. 1989]: 

JV 

^ h U^j (alaj + H.c. 



BHM 



i,j = i,i>j 
JV 

i=l 



(44) 



In this approximation of small anharmonicities, the total radial phonon number is conserved 
and coupling to axial modes is neglected. The third term in equation |44] has a contribution 
proportional to the phonon number squared which is responsible for the interaction. Its 
strength is scaled by [Deng et al. 2008] 

U = 2i-iyFTjlsw (45) 
The interaction depends on the strength and the relative position 6 of the standing wave and 
can even change its sign (see equations |42] and |45]l : it is maximum positive/repulsive for ions 
at the maxima of the standing wave {S — 0) and maximum negative/attractive for ions at the 
minima of the standing wave {6 = 1). This can be seen in figure[8]and|9l the small sinusoidal 
alteration from the standing wave makes the tip of the parabola slightly pointier for a minimum 
of the standing wave or chamfered for a maximum of the standing wave which results in a shift 
of the trapping frequency. The interesting consequence, however, are the differential shifts of 
the oscillator levels which are depicted in figure |9] the fourth order contribution from the 
extremum of the sinusoidal potential has always opposite sign of the quadratic contribution. 
Thus it is curved upwards for a maximum of the standing wave (meaning that the quadratic 
term is curved downwards) and shifts higher levels up, mimicking a repulsive interaction. A 
minimum of the standing wave has an upwards curved quadratic contribution; the fourth order 
is curved downwards, making the trap softer when walking away from the centre and reduces 
the energy of higher trap levels imitating an attractive interaction. 

The Bose-Hubbard model is well known in neutral atom physics [Fisher et al. 1989]. 
When the total phonon number A^ph is an integer multiple of the ion number A'ion the model 
predicts a transition from the superfluid state (dominated by tunnelling) to a Mott-insulator 
state (negligible tunnelling) with a constant phonon number per lattice site. This transition 
can be traversed by changing the ratio of tunnelling frequency and repulsive interaction 
energy (that is, changing trap frequency or standing wave intensity). For incommensurate 
phonon numbers that are no integer multiple of the ion number, a finite number of stages with 
constant phonon occupation can be anticipated, with the highest phonon populations at the 
bottom of the phonon confining potential, that is, at the trap centre. In neutral atom physics. 
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Figure 9. (colour online) Shifts of the harmonic oscillator levels due to the standing wave (not 
to scale). The centre image shows an undistorted harmonic trap with equally spaced energy 
levels. When a standing wave along x is added, such that a maximum is located at the trap 
centre (left hand side), the base level is lowered due to the quadratic part in the standing wave. 
In addition, the higher levels are raised due to the quartic part of the dipole potential. The 
levels spacing is no longer equidistant and mimics a repulsive interaction of the phonons. For 
an additional standing wave with a minimum at the trap centre (right hand side), the lowest 
level is raised due to the upward curvature of the dipole potential. The level spacing is reduced 
and also non equal and pretends a attractive interaction of the phonons. The grey dashed lines 
in the outer graphs again indicate unperturbed potential and energy levels. 



the superfluid-Mott insulator transition has been observed in beautiful experiments [Greiner 
etal. 2002]. 

The framework of Luttinger liquid theory [Luttinger 1963, Voit 1995] predicts algebraic 
decay of correlations between the number of phonons for the superfluid phase when phonon 
tunnelling beyond next-neighbours is neglected [Deng et al. 2008] 

= ("*"j> - (%) « K - jT' (46) 
and correlations that are nondiagonal in the phonon number decay as 

Ct:; = }, , / , (X \i - (47) 



with a oc y/UjtriQ given by the ratio of tunnelling and interaction at the central phonon 
density ng. In the Mott phase, however, both correlations decay exponentially with distance: 

with the correlation length ^. 

Phonon densities and fluctuations from numerical calculation using the density matrix 
renormalization group method (DMRG) [White 1993] were presented in [Deng et al. 2008] 
for linear traps and an array of microtraps. The calculations reproduced the superfluid Mott- 
insulator transition starting from the sides of the chain (due to larger ion separations, the 
tunnelling rate is lower at the ends). EspeciaUy for an array of microtraps the Mott phase 
features an almost flat mean phonon number per lattice site along the chain for commensurate 
ion and phonon number and the fluctuations are also flat and strongly reduced. For linear 
traps the inhomogeneous ion separations modulate phonon numbers and fluctuations which 
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are both more peaked at the chain centre. The correlations in the superfluid phase are found 
to be in agreement with ( 146 147 1 )) for short (linear trap) and intermediate (microtrap array) 
distances and signatures of finite size effects and inhomogeneous ion separations can be seen 
and are found to be more important in linear traps. In addition, long range correlations in 
Mott phase are found that can be attributed to the long range hopping terms in equation ( |40] |. 
For large repulsive interactions and incommensurate filling the system was found in [Deng 
et al. 2008] to form a Tonks-Girardeau-gas [Girardeau 1960, Lieb & Liniger 1963] which 
has been recently observed in neutral atoms [Paredes et al. 2004, Tolra et al. 2004]. Due 
to their repulsion, the phonons cannot get past each other and attain the hard core bosons 
limit and tunnelling is suppressed. The atoms/phonons behave neither purely bosonic nor 
fermionic, because they can occupy the same momentum state but not the same position in 
space. Mathematically, there is an exact one-to-one mapping between impenetrable bosons 
and non interacting fermions for a one-dimensional system. Numerical studies in [Deng 
et al. 2008] for an array of microtraps and half filling (A^ph = 1/2 A^ion) show a constant 
phonon number distribution together with algebraically decaying correlations. The exponent 
of the decay deviates from the value 1/2 which is expected for a Tonks gas with nearest 
neighbour interaction only. This is explained by mapping the BHM to a XY model with 
antiferromagnetic dipolar interactions. [Deng et al. 2008] also study attractive systems and 
use site dependent interactions to construct systems with highly degenerate ground states 
described by an XY model. 

4.2. Experiments 

For a string of 50 ions with a minimum axial separation of 5 pm and a radial trapping 
frequency v^^i about 70 times higher than the axial trapping frequency z/^ i the tunnelling 
rates tij were found to be comparable to the axial trapping frequency. For F = fiVx,i which 
is realistic with commercially available laser systems, the axial frequency is changed by a few 
ten percent and the system remains in the phonon number conserving regime and the phonon 
interaction energy U is approximately twice the tunnelling rate [Deng et al. 2008]. 

The interesting quantities are the number of atoms at each lattice site, corresponding 
to the number of transverse phonons at each ion, their fluctuations and their correlations. 
Phonon occupation can not be measured directly but can be mapped on the internal state and 
the ions' resonance fluorescence is recorded. The mapping is accomplished by sideband- 
transitions with carrier Rabi frequency VL - note that a sideband doesn't refer to the axial 
mode as usual but to the radial mode along the x axis here . Since the strength of the sideband 
transition depends on the expectation value of the phonon number, the transition matrix 
element, and therefore the vibrational state can be extracted from repetitive measurements 
[Meekhof et al. 1996fl, Meekhof et al. 1996^, Leibfried et al. 1996]. If the ion is not in a Fock 
state, each phonon population contributes to a sivl' {\/nVlt) term in the temporal evolution, 
which is accordingly given by 



with P{n) being the probability of having n phonons. Different frequency components in 
the temporal evolution can be extracted from Fourier analysis to yield the P„ [Meekhof 
et al. 1996fl, Meekhof et al. \996b, Leibfried et al. 1996] and reveal directly and transition 
from superfluid to Mott insulator state: For a Mott insulator, all ions are in the same vibrational 
state and therefore have the same transition matrix element (assuming the Number of phonons 
being an integer multiple of the ion number). Exciting sideband transitions on all ions 
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simultaneously brings all of them in the same internal state, so for appropriately chosen (tt- 
pulse) times the whole chain will switch from dark to bright. For a superfluid, however, 
the phonon expectation value is not a constant along the chain, and there is no universal tt- 
pulse and the ions will never undergo common oscillations from the dark to the bright state 
when varying the interaction time. If desired, the whole density matrix can be deduced from 
quantum state tomography which, of course, becomes tedious or, for all practical purposes, 
even impossible for long chains. 

4.2.1. Superfluid-Mott insulator transition and creation of a superfluid phonon state by 
adiabatic evolution This experiment is similar to the superfluid-Mott insulator transition by 
[Greiner et al. 2002] obtained by ramping the optical lattice depth. A clear demonstration, 
as proposed in [Porras & Cirac 2004^] would start with an ion chain cooled to a state with 
zero radial phonons, the "radial ground state". The optical lattice is ramped up to a high 
intensity U » Uj such that the ions sit on a node of the standing wave and the phonons 
experience a repulsive interaction. A defined number of phonons is introduced by sideband 
transitions, such that the total phonon number is an integer multiple of the ion number (for 
example by a tt pulse on the blue sideband). If now the phonon number is mapped back to 
the internal state by sideband transitions, a collective switching from the internal bright to the 
non scattering state should be observed, that is, the time for a tt pulse is identical throughout 
the whole chain). Alternatively, the interaction U is adiabatically reduced down to a final 
strength U f (the standing wave intensity is lowered), such that the system remains in the 
ground state. At a critical value Uf « tij the system undergoes a transition to a phonon 
superfluid. When measuring phonon occupation numbers by mapping them to internal states 
by sideband transitions, the phonon numbers can be expected to be inhomogeneous, due to 
the phonon trapping potential and non negligible tuimelling. 

4.2.2. Bose-Einstein-condensation by evaporative laser cooling Evaporative laser cooling 
for phonons in ion traps sounds hke an oxymoron because in neutral atom physics evaporative 
cooling terms a process free of laser scattering, where only the highest energy atoms can leave 
the trap [Davis, Mewes, Joffe, Andrews & Ketterle 1995] and take more than the average 
thermal energy with them which, provided sufficient fast thermalization) results in a cooler 
cloud. The analogy is, that in the ion trap with an axial phonon confining potential high 
phonon occupations (labelled high energy phonons in [Porras & Cirac 2004^?]) can more 
easily tunnel to the ends of the chain, analogous to hot atoms, and when laser cooling is 
applied to these ends exclusively, the phonons are removed or evaporated. 

A demonstration would start with a Doppler cooled ion chain with a given number 
of phonons per site. Applying laser cooling at the ends of the ion chain removes or 
evaporates high energy phonons (or phonon occupation numbers) from the top of the phonon 
confining potential. A small interaction U <^ tij allows for a phonon-phonon interaction and 
therefore for thermalization, equivalent to the need of a non vanishing scattering length for 
thermaUzation of neutral atom clouds and finally Bose-Einstein condensation of phonons can 
be observed. Since this is a ID system, BEC cannot take place in the thermodynamic limit 
A'^ — ^ oo at finite temperature, but for a given in number there will be a non-zero transition 
temperature. In contrast to the Mott insulator (phonon Fock state), the signature of this phase 
transition is more difficult to extract. One can deduce phonon occupation distribution from 
series of fluorescence images after sideband transitions of varying time. Evaporative laser 
cooling for phonons should reveal an abrupt transition from a thermal phonon occupation 
distribution to the BEC distribution, which corresponds to a distinction of a Gaussian and a 
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Thomas-Fermi momentum distribution for neutral atoms. 

It could be interesting to modify the axial trapping potential in microtraps such that 
the equilibrium positions are really equidistant and the phonon trapping potential, probably 
corresponding to unconfined phonons and a demonstration Bose-Einstein-condensation of 
free particles, similar to experiments by [Meyrath et al. 2005]. 



4.3. Frustrated XY models 

For large anharmonicities, strong repulsion, and low filling factor A'ph <C Mon the occupation 
by two phonons becomes unlikely and the system attains the hardcore Boson limit. The 
phonons can be mapped to = ^ spins via the Holstein-Primakoff transformation as 
(aj, 5'+, Oq — > and — > -S*^ + i). The resulting Hamiltonian is an XY model 
with dipolar interactions and a site dependent potential [Schmied et al. 2008]: 

Hs^2Y,to.,p (S"^ SI + SI Si) +Y.^o.Sl (50) 

{a,p) a 

Hamiltonians like theses are interesting because they allow to investigate frustrated XY 
models in one and two dimensions. Frustration is a phenomenon relevant in condensed 
matter physics and describes the situation where, due to the lattice geometry or competing 
interactions, the interaction energies cannot be simultaneously minimized which can result to 
highly non-degenerate ground states with non-zero entropy, even at zero temperature. Such 
frustrated XY models can now be investigated in ID and 2D system, that is, in linear chains 
and planar zigzag structures in linear ion traps and planar structures in surface traps. 

The ID case leads to the S* = ^ XY antiferromagnet with dipolar interactions, which 
becomes exactly solvable if only next neighbour interactions are considered. By weakening 
the transverse confinement along one direction the linear chain is no longer the energetically 
favoured pattern and the ion arrange themselves in a planar zigzag ladder which can be viewed 
as two linear chains relatively displaced to each other. By variation of the radial trapping 
potential, the properties of this two-dimensional structure, namely the zigzag amplitude and 
thus the ratio between the dominant inter (ti) and intra-chain (t2) coupling can be tuned. 
The zigzag amplitude ^, measured in units of the inter ion spacing, is limited by the fact 
that for a too weak transverse confinement the ion chain spontaneously flips into helical 
order (around ^ ~ 0.965). This in turn limits the accessible range of coupling ratios when 
phonons perpendicular to the zigzag plane are considered. But when exploiting phonons in 
the ladderplane any |i2/ii| > 1/8 ratio can be accessed. For low zigzag amplitudes and intra- 
chain couplings, the ladder is found to exhibit long range antiferromagnetic Neel order and 
undergoes a transition at ^ w 0.461 to spiral order which shows both spontaneous magnetic 
order and chiral order 

The chirality can be defined as 

= 4 (5f Sf+i - 5,f (51) 
and the chiral order parameter is the averaged chiral correlation over all ion-pairs in the ladder: 

= ^ _ ^_ («''^'+a) 

L-2 (52) 

A=-(L-2) 

When envisioning an experiment measuring this chiral order, we first have to remember 
that the spins in |5T| are not real spins but phonons mapped back via the Holstein-Primakoff 
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transformation and that the probability of having double phonon occupation is negligible. The 
expectation values for the phonon occupation are thus given by the probability of having single 
phonon occupation. To measure chirality and the chiral order parameter, we have to measure 
phonons in both transverse modes for each ion and their correlations by sideband transitions. 
When the measurement is repeated for different radial confinements and thus coupling ratios 
the transition from Neel to sprial order can be observed. Another interesting feature to look at 
is a yet not understood "reorientation transition" transition around ^ w 0.8 found in numerical 
calculations of up to 20 ions [Schmied et al. 2008] which gets more pronounced for larger 
ion numbers. Furthermore, the convergence of the numerical algorithm used in [Schmied 
et al. 2008] shows poor convergence for equal inter and intra chain couplings and thus a 
measurement of the phonon correlation function in this regime could be interesting. 

Even more variations are possible in extended planar lattices, as, for example triangular 
lattices as considered in [Schmied et al. 2008]. Rotating planar triangular Wigner crystals 
have been observed in Penning traps [Mitchell et al. 1998] and again they can be viewed 
as displaced linear chains such that it is possible to define a ratio of inter and intra chain 
coupling. A variety of phases has been found when the orientation of the a co-rotating 
standing wave, and therefore vibrational motion is changed relative to the lattice normal 
[Schmied et al. 2008]. Another approach to implement such lattices is by means of surface 
traps [Chiaverini & Lybarger 2008, Seidelin et al. 2006]. The equilibrium ion positions 
and therefore the ratio of inter and intra chain couplings is given by the effective potential 
generated by the trap. A variation of the ratio of relevant couplings could be investigated 
either in a series of surface traps, each designed for a particular coupling ratio, or by electrode 
geometries which allow for a anisotropic deformation of the triangular lattice. It should be 
pointed out that such planar traps would in principle allow to implement any desired lattice 
geometry. 

5. Axial standing wave: The Frenkel-Kontorova-Ion-Model 

The Frenkel-Kontorova-Model (FKM) is a model describing a hnear chain of particles with 
harmonic next neighbour interactions which is exposed to a sinusoidal potential [Frenkel 
& Kontorova 1938]. The system of interest might remind the reader of the Hubbard- or 
Bose-Hubbard-model as discussed in section 14.11 however, the Frenkel-Kontorova-Model 
is simpler, being classical and one-dimensional, and, on the other hand, focusses more on 
classical nonlinear dynamics and chaos, and questions of integrability and the persistence 
of quasi-periodic motion [Braun & Kivshar 1998]. Despite this simplicity the FKM is 
able to describe a variety of physical phenomena as crystal dislocations, commensurate- 
incommensurate phase transitions, epitaxial monolayers on a crystal surface, magnetic chains, 
and fluxon dynamics in Josephson junctions (see [Braun <& Kivshar 2004], and references 
therein). The FKM is one of the first examples in solid state physics where a one dimensional 
chain can be used to model an extended two dimensional defect in a bulk [Braun & 
Kivshar 1998]. 

5.7. Formalism and numerical calculations 

In the basic model the interactions are harmonic and restricted to next-neighbours and the 
external potential is sinusoidal (see figure [TOli. Without this external potential the dynamics 
of the chain can be understood in terms of normal modes and is quasi-periodic. When the 
potential is switched on but the perturbation is small, that is, below the critical perturbation 
strength, the system is in the so called sliding phase and the chain can oscillate in the lattice 
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a) 




b) 




c) 




Figure 10. a) basic Frenkel Kontorova model: a one dimensional chain of atoms with next 
neighbour interactions is exposed to a sinusoidal potential; figures b) and c) show possible 
spatial configurations: the upper configuration b) is unstable - the total energy, which is 
a function of the coordinates of all particles has a saddle point, whereas configuration c), 
obtained only by a minuscule sideways shift is stable, con'esponding to a minimum of the total 
energy (after [Braun & Kivshar 1998]) 



and the ion positions follow a Kolmogorov-Arnol'd-Moser curve [Braun & Kivshar 2004]. 
For higher perturbations, the system undergoes an Aubry analyticity breaking transition 
[Aubry & Le Daeron 1983] and the ion chain starts to be pinned by the lattice and the positions 
form a devils staircase corresponding to a fractal Cantor set [Garcia-Mata et al. 2007]. The 
ground state is unique for all perturbations, but in the pinned phase, a large number of states 
exists very close to the ground state like in a fractal spin glass [Zhirov et al. 2002]. 

Various extensions and modifications of the basic Frenkel-Kontorova model have been 
investigated, for example arbitrary on site potentials and anharmonic interactions between the 
particles [Braun & Kivshar 1998]. In the following we will focus on the Frenkel-Kontorova- 
lon-Model (FKIM) which is a classical picture for a laser-cooled chain of ions trapped in a 
linear trap exposed to a standing wave [Garcia-Mata et al. 2007]. For a sample of ions in a 
linear trap with a much stronger radial than axial confinement, the ions arrange themselves 
in a linear chain. When the radial motion is frozen due to laser cooling, the chain can be 
described in good approximation with a one dimensional Hamiltonian. Then the interaction 
of the particles is given by the long range Coulomb force and the external potential is a sum 
of the effective sinusoidal potential from the standing wave and the effective harmonic axial 
trapping potential. The Hamiltonian is given by 



H 



N 

E 




K cos{2'KZi/d) 



(53) 



47reo \zi 



The Hamiltonian above corresponds therefore to the general ion trap Hamiltonian given in 
equation ( [36l ) where only the part in z-direction is considered and an sinusoidal potential 
of strength K and periodicity d is added. A basic FKM Hamiltonian can be obtained 
by neglecting restricting the interaction sum to next neighbours and expanding the sum of 
interaction and trapping potential in the deviations from the equilibrium positions up to second 
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Figure 11. Frenkel-Kontorova ion model: The potential curvature, the sine ripple amplitude 
and lattice constant are plotted for values corresponding to the Auhi-y analyticity breaking 
transition at golden mean density for a chain of 12 ions with a centre separation of 5 ]im. A 
CCD image of a string of 12 ^^-^Yb+ions is scaled to match the conditions for the potential 
plot and indicates the unperturbed ion positions. 



order. In [Garcfa-Mata et al. 2007] the authors investigate the dimensionless form of the 
FKIM-Hamiltonian above which reads 



The ion positions z are given in units of the reduced lattice constant d — d/2TT and 
the energy E as well as the strength of the sinusoidal potential K are measured in units of 
e^/47reo(i. The modified angular frequency is obtained as = ATi€()rav1 ^d^ /e^. In the 
quantum case the momenta are given by Pi = —ificsdxi with an effective Planck constant 
Ues given as 



which depends on the lattice constant which allows to scale and investigate quantum effects 
as we will see later. 

First the authors compare numerical solutions to the classical ground state problem in the 
FKIM to solutions of the FKM with only next neighbour interactions. The numerical studies 
are performed for up to a few hundreds of ions at the golden mean density (at the center of the 
ion chain), where the average occupation of a lattice site is (\/5+l) /2 and results in the golden 
Kolmogorov-Arnol'd-Moser (KAM) curve for the ion positions commonly used for studies of 
the Aubry transition [Aubry & Le Daeron 1983, Braun & Kivshar 2004]. The FKM leads to 
dynamical recursive maps for the calculation of equilibrium ion positions which are found to 
describe the centre 1/3 part of the ion chain appropriately. By comparison of the solutions 




(54) 




(55) 



Quantum Simulations with Cold Trapped Ions 



34 



for FKM and FKIM equilibrium positions, the influence of long range interactions and the 
trapping potential can be mapped out and is found to be negligible in the centre region. Here, 
the ion density is highest and almost constant and allows for screening of the trap potential 
by nearby ions, similar to the Debye radius in plasma physics. In the outer regions visible 
deviations in the equilibrium positions are found, as the ion density is lower at the ends of the 
chain. 

As for the FKM, for small depths of the sinusoidal potential (K < 0.05), the chain is in 
a sliding phase. The hull function, which compares the unperturbed equilibrium positions to 
the perturbed ones modulo 2tt is continuous and the phonon mode spectrum has a sound like 
form. At a critical strength of the standing wave {Kc ~ 0.05), the system undergoes an Aubry 
analyticity breaking transition [Aubry & Le Daeron 1983]. Above the critical lattice strength, 
the ion chain is in a pinned phase, their positions are described by a devils staircase, which 
is a fractal cantor set [Peitgen et al. 2004], and the phonon spectrum shows a gap. Another 
similarity to the FKM is that in the pinned phase, the FKIM has properties of a spin glass, and 
features an enormous number of stable equilibrium configurations energetically close to the 
ground state [Zhirov et al. 2002], despite the ground state being unique. 

So far, all considerations have been purely classical. Quantum effects were in 
investigated in [Borgonovi et al. 1989] and the quasi-degenerate configurations close to the 
ground state are found to become important and tunnelling between them leads to non- 
trivial instanton excitations [Garcia-Mata et al. 2007]. The quantum case is studied in 
[Garcfa-Mata et al. 2007] using the Quantum-Monte-Carlo (QMC) approach [Grotendorst 
et al. 2002, Nightingale & Umrigar 1998]. To see this structural changes, the authors look 
at the power spectrum of the position distribution (in the classical case, the z„ are the ion 
positions): 



From these numerical solutions a phase transitions from a pinned instanton glass to a 
sliding phonon regime is observed, depending on the value of the effective Planck constant 
hcs [Garcfa-Mata et al. 2007] which can be varied by changing the lattice constant and thus 
the influence of tunnelling, and the transition taking place at TicS ~ 1- For low values of 
hes the formfactor shows discrete resonances which are equidistant in the sliding phase, and 
and integer k showing that the ions are pinned by the optical lattice. The pinned phase also 
shows a quantum transition when hcs is increased: for low values hcff, the pinned phase is 
not destroyed and discrete resonances are still found. For larger values of h the formfactor 
bacomes continuous and the transition occurs at a value TicS ~ 1- 

5.2. FKIM Experiments 

Experiments for investigating Frenkel-Kontorova-physics require a sinusoidal potential along 
a one-dimensional ion chain. This can be implemented by a far off -resonant standing wave 
along the weak trap axis as indicated in figure[T2] Being far off resonant means that changes of 
the internal state due to this lasers is unlikely within the duration of the xperiment. Following 
[Garcfa-Mata et al. 2007] the ion density at the centre should be the golden mean density, that 
is, the lattice periodicity should be by a factor (\/5 + l)/2 larger than the ion separation at the 
middle of the chain. This can easily be obtained even at a fixed wavelength changing the angle 
between the intersecting beams that form the dipole potential or, alternatively, by changing 
the trap frequency. To obtain the golden mean density for ions separated 5 pm, the lattice 
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a) b) 




Figure 12. Schematic (not to scale) of the experimental setup suitable for simulating Frenkel 
Kontorova physics and observing an Aubry transition. The standing wave is applied in the 
axial direction, either by optical access through ring shaped endcap electrodes, or by crossing 
two beams at an angle, such that the counter-propagating component is aUgned parallel to the 
weak trap axis. Part a) shows the view along the ion chain, part b) shows the setup as seen by 
an imaging detector used for measuring the ion fluorescence and for the determination of ion 
positions and formfactor 

periodicity is approximately 3.1 |jm corresponding for example to two beams at A = 1064 nm 
each enclosing an angle of a « 80° with the weak trap axis. 

The position of the lattice has to be fixed relative to the ions which means that the 
relative phase of the two counter-propagating running waves has to be constant. This requires 
stabiUty and potentially active feedback on mirrors. To maintain a constant depth of the 
sinusoidal potential along the chain, the Rayleigh range or would ideally be large compared 
to the chain length, meaning that the waist would be large. The required strength of the 
standing wave Kc to observe an Aubry transitions corresponding to a depth of 0.6 K ions 
separated by 5 [im which is much deeper than in neutral atom experiments. On the other 
hands the trap frequencies are much higher in ion trap experiments and thus the time scales 
for a measurement are shorter, allowing for lasers, that are closer to the atomic resonance and 
in turn requiring less power. The challenging power requirements can be lessened, when the 
light cycles in a resonant cavity, however, with a more demanding stability. 

The effective Planck constant as given in equation ( fSST l contains the lattice constant and 
can in principle be changed by adjusting either the wavelength of the light that forms the 
sinusoidal dipole potential or the angle between the intersecting beams. To maintain the 
golden mean density, the trap frequency would be varied simultaneously. For parameters 
given above, however, the effective Planck constant is approximately 3 • 10—5 (calculated for 
Ca ions [Garcia-Mata et al. 2007]) and due to the weak dependence on the lattice constant 
(see equation ISSTi ion trap experiments are likely to be restricted to the semi-classical regime 
with ?ioff ^ 1, hardly accessible for QMC calculations [Garcfa-Mata et al. 2007]. 

The experiment would start with an ion chain cooled to the ground state. Ramping 
the standing wave along the weak trap axis adiabatically up to the desired final strength 
ensures that the system remains in its ground state and the equilibrium positions and thus 
the formfactor can be measured. 

Beside of the general interest in Frenkel-Kontorova physics, an axial standing wave as 
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discussed above might also prove to be useful in the context of quantum information with ion 
traps since the strong gap in the phonon spectrum found in the pinned phase might allow for 
the protection of quantum gates against decoherence [Garcia-Mata et al. 2007]. 

6. Quantum Fields and relativistic effects 

6.1. Introduction 

The way to prove a physical theory is to compare its predictions with the observations. 
This has been the case with the cosmological models where effects like the red shift or the 
homogeneous background radiation constitute observations in accordance with a description 
of an expanding universe and the particle creation. But to perform a truly cosmological 
experiment is out of the reach of our current technology, in particular one that reproduces 
the early stages of the universe characterized with an extremely high energy density and a 
rapid expansion. 

In other cases like with relativistic particles, it is possible to make experiments were 
electrons for example are driven to almost the speed of hght. But some interesting effects are 
still out of reach of our measurement capabilities. For example the Zitterbewegung, predicted 
from the Dirac equation. It states that a freely moving spin- 1/2 particle, in absence of 
external potentials, is subject to a helicoidal motion around the main direction of propagation. 
Another counterintuitive effect, the Klein paradox [Klein 1929] states that an electron could 
transmit unimpeded through a potential barrier. But this effect has not been observed yet with 
elementary particles. 

6.2. Simulating the Unruh effect 

6.2.1. Methods Following the work of Unruh [Unruh 1 976] it was shown that in the vacuum 
of a Minkowski (flat) universe an accelerated observer would measure a thermal spectrum 
of particles whose temperature is proportional to his acceleration. Later this effect was 
generalized by Gibbons and Hawking [Gibbons & Hawking 1977] for a curved de Sitter 
space. The de Sitter space is a massless curved spacetime that is isotropic and homogeneous 
in space, and which looks like the Minkowski space locally [Bergstrom & Goobar 2004]. The 
de Sitter space is also a cosmological model that describes an exponential expansion of the 
universe due to a non vanishing and positive cosmological constant A. 

The scalar field describing a massless inflating universe with a single global mode can be 
imitated by the quantized motion of a single ion around its equiUbrium position. The phonon 
excitations of the ion would simulate the field quanta of the universe. The ion itself can be 
used as a detector by coupling its vibrational state to an easily readable internal state. 

In the case of an expanding de Sitter universe a comoving quanta field detector would 
experience a Doppler frequency shift of the form iy{t)^ = cxp (— VSAf) [Birrell & 
Davies 1984]. The spectrum measured with this detector is equivalent to the one measured 
by a thermally excited inertial detector with a temperature T = %\/^K/kB [Gibbons & 
Hawking 1977]. 

In the case of a single ion trapped in a time dependent potential that varies in the 
form v{t) = i/Qexp (— Kt), a similar effect can be measured. A pseudo-temperature of 
ksT = hK/2-K is associated with the average occupation number of the vibrational energy 
levels. 

The coupling between internal and vibrational states of the ion can be created by a 
laser field described in the dipole and rotating wave approximation by the interaction-picture 
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Hamiltonian [Alsing et al. 2005, James 1998], 

Hi = hnkz{t) [a-e-'"^^ - a+e'^^] . (57) 

Here, and <t^ are the raising and lowering spin- 1/2 operators, A is the detuning of the 
laser relative to the atomic transition and fl is the Rabi frequency, respectively. The laser 
wave vector k is pointing along the direction in which the ion's vibrational motion is excited 
(here the direction). We assume an anisotropic trap (e.g., a linear trap) where, by choosing 
the detuning A and the wave vector k, the experimenter can control which vibrational motion 
to excite. 

The evolution of the position operator is deduced by solving the equations of motion 
of the ion in its harmonic trap neglecting its coupling with the Ught field. In the Heisenberg 
representation and with 

H=f + ^.{trz' , (58) 
2m 2 

we have [Menicucci & Milburn 2007] 

*rff 1 
z = -t[H,z] = — 

, ^ (59) 

Pz = = mv{t)'^z 

After the opening of the harmonic potential v{t) during a time T in which the interaction 
is on, a vibrational phonon can be mapped into the excited 1 1) internal state. The probability 
of exciting the internal state of the ion calculated with first order perturbation theory and in 
the slow K <C i^o and long-time t'oe^"^ <C k, frequency variation for the trap potential is 
[Menicucci & Milburn 2007], 

P(1)(A) = ^ dt, dt, (0| Hiih) |1) (1| Hj{t2) |0) 

When the trap is opened to simulate the expansion of the universe, the ion is used as a 
detector to measure the thermal spectrum from the phonon distribution. A particular case is 
when the laser is in resonance with the red sideband (A w ly). The ion acts like a phonon 
detector that absorbs the measured particle (together with a photon) in order to produce an 
internal excitation of the ion. When the laser is in resonance with blue sideband (A w — i^), 
the ion acts like an unconventional detector that emits a phonon to get excited. The ratio of 
the excitation probabilities gives the experimental signature of the ion thermal distribution 
[Wineland & Itano 1979, Monroe, Meekhof, King, Jefferts, Itano, Wineland & Gould 1995], 

where h is the average phonon occupation value. From the excitation probabilities of equation 
|60]the ratio between the red and the blue sidebands is given by [Menicucci & Milburn 2007], 



= e-^'^'^/'^. (62) 



P(i)(-A) 

Therefore the temperature associated with the exponential opening of the trap potential is 
given by T = fiK/2'iTkB- This temperature is equivalent to the Gibbons-Hawking temperature 
T ^ hVSA/kB [Gibbons & Hawking 1977]. 

The rise of the temperature of the trapped ion is due to the trap opening. At the beginning 
of the trap opening the evolution is adiabatic since the trap frequency obeys the condition 
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<C v. Indeed, the expansion constant was selected small compared to the initial trap 
frequency k <C i^o- At a later time T chosen following the condition i^^e^'^'^ ^ k the 
evolution turns to be non adiabatic. Indeed, the trap frequency is going to be so weak that even 
a small rate of change would be non-adiabatic. The ion that initially was in the ground state 
is excited to higher vibrational levels because of the non-adiabatic reduction of the confining 
potential. 

6.2.2. Experiment In this case we measure the temperature associated with the ion phonon 
distribution. This temperature is related to the Gibbons-Hawking temperature. 

To perform this experiment a single trapped ion is cooled to the vibrational ground state 
level. Next the ion is pumped into its electronic ground level that is related to the vacuum state 
of the universe. Then, the trap potential, initially at the frequency vq, is exponentially reduced 
v{t) = exp {—Kt). During the opening of the trap, Raman laser beams between hyperfine 
levels or a narrow frequency laser tuned to a quadrupole transition are set in resonance with 
the red or the blue sideband. It should be noticed that the couphng field must be in resonance 
with the sideband during the trap opening. After some time, the ion internal state is measured 
by fluorescence using the shelving technique with the help of an additional laser beam in 
resonance with the excited internal state. 

The trap potential is then raised and the sequence can be repeated for the opposite 
sideband. The height of the red and blue sidebands are compared and the temperature 
dependence on the trap expansion coefficient k is measured. 

A typical value of the initial trap frequency vq ~ IMHz can be used to perform this 
experiment. Therefore the expansion coefficient can be selected like k ~ IkHz and the 
expansion time T of several ms. These are typical parameters in current ion trap experiments. 

Difficulties that may arise when performing such an experiment are associated with 
unwanted heating of the ion caused, first, by a possible unintentional displacement of the 
trap minimum and by anharmonic terms in the trap potential, and second, by fluctuating 
patch potentials. The former could be solved by compensating the anharmonicities and the 
displacement of the trap potential by tuning appropriate compensating fields for every value of 
the time- varying trap frequency. The latter could be overcame by performing the experiment 
during a short enough time so that the heating can be neglected. 

6.3. Simulating the cosmological particle creation 

6.3.1. Methods The standard model describes the particle creation from the vacuum 
fluctuations during the early stages of the universe. Since the expansion of the universe is non 
adiabatic, an initial ground state will be excited evolving into a squeezed state of entangled 
pairs of particles. The created particle pairs are simulated by entangled phonon pairs created 
by the excitation of the ground ionic motional level as the trap potential is non-adiabatically 
varied. 

One can use a real massless scalar field (6 as a an approximation to describe an expanding 
or contracting universe. Using the Friedman-Lemaitre-Robertson- Walker metric [Birrell & 
Davies 1984], ds^ = a^{t)dt^ — a'^{t)dr'^ one finds that the wave equation describing the 
universe reads [Schiitzhold et al. 2007], 



where the Ricci (curvature) scalar 3? is coupled to the field via the dimensionless parameter 
(. Each mode k represents a harmonic oscillator with a time dependent potential a**(f)k^ + 




(63) 
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In the case of a chain of ions trapped in a transversally strongly confining Unear trap with 
a time dependent axial trap frequency (t), the position Zj of the /th ion follows the equation 
of motion 

■• , 2/,^ '^sign(i-j) 

z, + V, {t)zi =7 2^ 7 72 (64) 

This equation can be solved classically by the use of the scaUng ansatz, 

z,{t) = (65) 
where z° are the initial equilibrium positions. The scaling factor h{t) evolves following 

S + (06) 

The quantum fluctuations are introduced by splitting the position operator Zi into its classical 
evolution h{t)zf and a quantum fluctuation term 6zi{t), 

Ut)=b{t)z^ + dz,it). (67) 

Replacing equation|67]in equationlMland linearising later over 6zi (t) (since the displacement 
of the ions can be assumed small) one obtains after normal mode expansion [Schiitzhold 
et al. 2007], 

^ + [-fW + ^])^^. = (68) 

for the phonon modes k. One can recognize the similarities between equation|63]and equation 
|68] The scalar field 0k behaves in the same way as the quantum fluctuation of the ions 
(phonon) (Jz^ and each mode k of the field has a direct analogy with the phonon frequency of 
the ions i/^. 

It has been shown that a non-adiabatic change in the trap potential can create squeezed 
states [Heinzen & Wineland 1990]. 



6.3.2. Experiment The simulation can be performed in a system composed of, for example, 
two hyperfine states |0) and |1) of an ion. The states should be coupled by a two-photon 
Raman transition or by an RF field in the presence of a magnetic gradient. Notice that the pairs 
of Raman beams or the RF field not only couple the two hyperfine states |0) and 1 1) but also 
the vibrational levels of the ion \n) — {\n = 1) , |n = 2) , ...). For this simulation one uses 
three different coupling fields, one coupled to the second red sideband (a) |0, ?i + 2) |1, 
another to the first red sideband (b) |0, 7i + 1) 1 1, n) and the last one to the carrier transition 
(c) |0, n) <-> \l,n). An additional level |e) is used for fluorescence detection of the population 
in level |0). The laser pulse (d) |0) |e) drives a dipole transition that does not resolve the 
vibrational levels as the Raman beams or RF field do. 

A non-adiabatic increase of the strength of trap potential confining the ion will simulate 
the universe expansion. After non-adiabatically changing the trap potential, the ions motional 
state will be a squeezed state where only even (i.e., symmetric with respect to the trap 
minimum) vibrational levels contribute. Thus, the signature of the squeezed state is the 
population of even vibrational levels. Measuring the population of the \n = 2) vibrational 
level should be sufficient, since the population in higher vibrational levels is expected to be 
negligible. In order to measure the population of the n = 2 state one performs a series of 
pulses described in figure[T3]that couple the motional state to the easily readable internal state 
of the ion. 
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Figure 13. A possible energy level diagram and the transitions used to measure the population 
in the n = 2 level. If the ion is found in level |0, n = 2) after the trap contraction, then 
applying a tt second red sideband pulse (c) followed by a tt carrier pulse (a) brings the ion to 
the state |0, n = 2). When the ion is in this latter level, fluorescence can be seen by using the 
detection laser {d). On the contrary if the ion is initially in |0, n = 0) or |0, n = 1) the first 
pulse (c) would not change the state of the ion since it does not have enough energy to be in 
resonance with the |1) energy level but the subsequent carrier pulse (c) transfer the population 
to the corresponding vibrational levels of Then no fluorescence can be seen when the 
detection laser (d) is applied. Note that the transitions a and c should be narrow enough not to 
couple the adjacent vibrational levels and that the detection laser drives a dipole transition that 
couples many of the vibrational levels. 



Since cooling is more efficient for strong confinement, it is advantageous to begin with 
a strongly confining potential and cool to the ground state of motion. Then, after an adiabatic 
lowering of the trap potential it is ramped back up non-adiabatically to its initial strength. 

To discriminate against any classical heating, i.e. parametric or other classical 
perturbations, the probability of populating the n — 1 level can be measured in a similar 
way as shown in figure [13] The comparison between both probabilities gives the proof that 
a squeezed state is generated, if a higher probability for populating the n — 2 level than the 
n = 1 is obtained. 

6.4. Simulating the Dirac Equation 

6.4.1. Methods Relativistic effects related to the solution of the Dirac equation can be 
simulated in trapped ions systems as proposed in [Bermudez et al. 2007] and [Lamata 
et al. 2007], in particular phenomena like the Zitterbewegung (helicoidal motion of the free 
Dirac particle as a consequence of the non-commutativity of its velocity operator components) 
or the Klein paradox (transmission of a relativistic particle of mass m through a potential edge 
of height V > 2mc^ with nearly no reflection). 

In order to simulate the Dirac equation in 3 + 1 dimensions for a spin- 1/2 particle in a 
single trapped ion, L. Lamata et al. proposed to design a Hamiltonian that reproduces the same 
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dynamics. For a free electron in three dimensions the Dirac formalism uses a bispinor that is 
a vector valued wave function of four components. In order to represent all the components 
one can use a four level system with states \a), \b), \c), \d) that correspond to four internal 
levels of the ion. The bispinor is defined as [Lamata et al. 2007], 

\^) := *a \a) + *h \b) + ^-e |C) + \d) - 



(69) 



The dynamic of the system is generated by the pairwise coupling of the internal levels with 
the center of mass vibrational levels whose energies are determined by the trap frequencies 
Vx, Vy, Vz. Again, three common types of interaction are used. A carrier interaction 
coherently couples two ionic internal levels without changing the external motion of the 
ion for example \a.n) ^ |c, «-). A Jaynes-Cummings (JC) and an anti-Jaynes-Cummings 
(AJC) interactions resonantly couple the ionic internal with the vibrational levels. The former 
induce the red sideband transition and the later the blue one for example |a, rt + 1) ^ |c, n) 
and |a, n) ^ |c, n + 1) respectively. The corresponding interaction Hamiltonians under the 
rotating wave and Lamb-Dicke approximations are. 



(70) 



Here again and are the raising and lowering ionic spin- 1/2 operators, and a 
are the creation and annihilation operators associated with the motional state of the ion and 
Tj — yJTiW/2Mv is the Lamb-Dicke parameter The Rabi oscillation frequencies for the blue 
and red sideband transitions were set to the same value Cl and to Q, for the carrier transition. 



When only the carrier coupling is applied, with the appropriate choice of the phase 4> one 
obtains, 

H^^ = hQ.jaj (71) 

with j = {x, y} and aj being the Pauli matrices. 

When the blue and red sideband couplings are used simultaneously, the resulting 
Hamiltonian is the sum of Hr and Hb from equation |70] By choosing the phases of the 
sideband fields like (pi, — (f)r — tt one obtains, 

= ±ihr]Anaj{a - at). (72) 

Here, A = ^Ti/iMv is the spread in position of the ground wave function. The operator 
(jj is the Pauli matrix along the j = {x,y} component of the spin. It is not related to the 
vibrational directions of motion but to the ionic internal transitions. 

The selection of the Pauli matrix ctx — cr^ + a~ or (jy — i{a^ ^ cr^) well as 
the sign of the Hamiltonian is performed by an additional condition of the applied red and 
blue sidebands phases. Even if a condition on their difference was already used, one needs a 
supplementary condition to completely define the phases. The sum of the phases of the applied 
fields = {(i>b + '?i'r)/2 determines the sign of the Hamiltonian, negative for 0^ = {tt, 3/27r} 
and positive for (j)^ — {0, I/Stt} whereas the Pauli matrix aj is selected along j = x for 
(f>s = {0, tt} and along j — y for (ps — {7r/2, 3/2tt} (see figure [T4]|. The phase (j>s itself does 
not modify the dynamics of the system when the previous Hamiltonian is considered alone. 
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but it is essential to build the Dirac like Hamiltonian composed of many elements couplings as 
in equation|72] The dynamics of equation|72]are based on a spin dependent coupling between 
the internal and vibrational states that can also be used for the creation of entangled states 
of spin and motion. This coupling has already been applied in previous experiments [Haljan 
et al. 2005, Sackett et al. 2000]. 




Figure 14. The selection of the phases of the red and blue sidebands (fir and t/ij, fields are 
important since they shape the form of Hamiltonian. In particular the sum of the phases 
03 = (0r + <^b)/2 defines the direction of the spin operator and also the global sign of 
the Hamiltonian l72l 

The Hamiltonian of equation |72] constitutes the basic block to construct the Dirac 
Hamiltonian in (3+1) dimensions. Since the spin-1/2 particle to simulate has a momentum 
with components in aU three spatial directions, one needs to consider Hamiltonian terms along 
all of these directions. The operators and a can be associated with the three normal trap 
frequencies and therefore with the motion along the three trap axes. One needs to distinguish 
a\ and for every axis £ — x,y, z. The same consideration should be made for the Lamb- 
Dicke parameter ry^, the spread in position of the ground state wave function and the 
Rabi frequency fig that also need to be defined along the direction of the applied interaction. 
The difference between the creation and annihilation operators can be written in terms of 
the momentum operator pi = ih,{al — ai)/ Ag. Following the previous considerations, the 
Hamiltionian of equation |72] can be rewritten as, 

HP^^ = 2rjeA,n,ajPi. (73) 

In order to simulate the Dirac dynamics employing four internal states of an ion, the 
transitions \a) ^ |c), \a) ^ \d), \b) <-> |c) and \b) ^ \d) are employed. The corresponding 
spin-1/2 operators of equation |73] need to be defined for the specific transition a^'^, a^'^, 
and a^^ and therefore the corresponding Hamiltonian from equation |73] can be written as 
IP'^ , E^"^ I ,,, H-P'^ , and ... In the case of the carrier interaction of equation [TTI 

aj{acy ajyady <7j{bc) <7j[bd) ^ ' ' 

two internal transitions are used, giving the following Hamiltonians i?crf(ac) H^^ijj^y 

With the correct selection of the direction and phase for the different fields, one can 
tailor the desired Hamiltonian as a linear combination of Ha^, and if^* yielding, for instance 
[Lamata et al. 2007], 



^ a^iad) <Tx(dc) ay(ad) <^v(bc) (74) 



Quantum Simulations with Cold Trapped Ions 



43 



It should be noticed that all the Rabi oscillation frequencies for the sideband and carrier 
transitions f2 and fi, the spread in position A and the Lamb-Dicke parameter 77 are selected 
equal in all the components of the previous Hamiltonian for reasons that will be explained 
later. This means that the fields employed to drive the various couplings should have the same 
relative strength for all the directions of space and the ion trap should be spherical. Therefore, 
one may obtain a Hamiltonian Hd that is symmetric with respect to all directions of space. 
The Hamiltonian Hd can be written in a 4 x 4 matrix representation in a base defined by the 
four internal states of the ion (see eq. |69] |, 

^ ^ V ir^mia ■ p) + inn 

where each element represents a 2 x 2 matrix. This Hamiltonian can be compared with the 
"supersymmetric" representation of the Dirac Hamiltonian [Thaller 1992], 

iin={ ' . . (76) 

\^ c((T • p) + ivnc j 

The speed of fight c and tfie electron rest energy m(? are related to tfie coupfing strengtfis 
and following, 

c = 277Afj and ttk? = hVL. {11) 

Tfie previous selection of symmetric couplings applied along all tfie directions of space 
represents in a simulation tfie isotropy of tfie velocity of ligfit. 

One can recognize tfie similarities between tfiese two Hamiltonians, making possible 
tfie simulation of a quantum relativistic evolution of a spin- 1/2 particle witfi an ion-trap 
experiment. 

If tfie dynamics created by tfie Hamiltonian of equation |75] were reproduced in an 
experiment witfi a four level system using Raman beams, it would require 14 pairs of Raman 
lasers, two for eacfi element of tfie form i/^j (see equationl73]l and one for eacfi carrier element 
H^i, (see equation ItTT). Of course, some simplifications can be applied, like modulating tfie 
frequency of tfie beams to generate tfie blue and tfie red sideband transitions witfi a single pair 
of beams. But it sfiould be stressed tfie need of controlling tfieir pfiases independently and tfie 
minimal requirement of fiaving at least two pairs of beams (modulated in frequency) witfi a 
resulting k-vector along eacfi one of tfie trap axes (see figure [TSl). 

In order to render tfie simulation acfiievable witfi current experiments, one can simplify 
tfie 3h-1 dimensional Hamiltonian of equation |76] and work witfi 2h-1 dimensions [Bermudez 
et al. 2007], or simplify even more to 1h-1 dimensions tfiat still reproduce tfie most important 
effects of tfie Dirac equation. Tfie resulting 1h-1 dimensional Hamiltonian becomes, after a 
7r/2 rotation around tfie x— axis tfiat cfianges (jy into [Lamata et al. 2007], 

H^j^'^ = 2T]An(T^p^ + hna^,. (78) 

Tfiis Hamiltonian operates on tfie Dirac "spinor" j^f^^^) — '^'^a^ \a) + \b) build witfi 
a positive (^'i^'') and a negative (^E'l,^'') kinetic energy component [Tfialler 1992]. Tfie term 
2riAil(TxPx results from tfie application of blue and red sideband fields. It is just tfie sum of 
Hi, and Hr witfi pfiases (pi, ~ ~7t/2 and <j)r — Tr/2. Tfie last term Tiflaz can be obtained, for 
instance, by applying an off-resonant laser pulse tfiat induces an ac-Stark sfiift. Because of 
tfie reduction to one spatial dimension, tfie trapping potential sfiould also be one-dimensional, 
tfiat is, a linear trap witfi strong radial confinement compared to tfie axial one would be well 
suited. 

Tfie experimental requirements for simulating tfie 1 + 1 Dirac Hamiltonian are 
substantially reduced and acfiievable in current experiments. 
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Figure 15. Possible laser configuration to build the Dirac Hamiltonian with Raman ti'ansitions. 
In order to build a symmetric Hamiltonian as the one of equation 1751 the wave vectors must 
be the same along the three trap axis. For simplification of the lasers requirements only one 
wavelength can be used. Then the four levels to perform the simulation need to be close 
enough in energy i.e. two Zeeman split Hyperfine levels as proposed in [Lamata et al. 2007]. 
Therefore the angles subtended by every pair of Raman beams must be the same. Additionally 
the beams need to be modulated with the correct frequency to drive the desired transitions. 



6.4.2. Experiment for the simulation of the Zitterbewegung By numerically solving the 
Dirac equation, it is possible to visualize peculiar behaviour of the probability density 
function. One of this peculiarities is the Zitterbewegung where the position of a particle 
oscillates even in the absence of external potentials. Usually this relativistic effect is not seen 
in experiments since the frequency is on the order lO^^Hz and the amplitude of order lO^'^A 
for a real particle such as an electron. 

When simulating the Zitterbewegung with trapped ions the frequency and amplitude of 
the oscillatory motion for an ion with average momentum po ^"e given by [Lamata et al. 2007] 



From equation |79] one can obtain a measurable output with frequencies ranging from zero 
to some Megahertz and amplitudes that can go to a thousand Angstroms varying the typical 
parameters of Rabi frequencies and trap depths. It can be noticed that the amplitude of the 
Zitterwebegung can be made independent of the Rabi frequencies if 77^! w indicating that 
this effect could be seen even with small laser fields. 

In order to observe the Zitterbewegung one needs to determine the momentum and 
the position of the ion for different times. In previous experiments [Meekhof et al. 1996a, 
Meekhof et al. 1996^] the state of motion of the ion was determined by measuring the 
population distribution of the Fock states. From the red sideband excitation probability 
as a function of time one can fit the different frequency components associated with the 




(79) 



rjh^nnA 



(80) 
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contribution of every vibrational level. Therefore, the determination of the distribution over 
Fock vibrational states requires many data points. Another possibility could be to measure 
the Wigner function by instantaneous measurements, even in the presence of a thermal bath, 
based on the proposal in [Santos et al. 2007]. 

To reduce the number of data points to acquire, one can just measure the expected 
value of the generalized quadrature = (ae^**^ — a^e"^)/2 as proposed by [Lougovski 
et al. 2006]. Indeed one does not need to measure the whole wave function to know the 
position or momentum of the ion, it is enough to measure the expected value. The angle 
(f) is related to the internal spin state I+0) = (|a) + e**^ \b)) where the ion is initialized. 
Additionally, a JC interaction needs to be applied (see Hr in equation iTOli which couples the 
vibrational state to the internal state where the measurement is performed. 

The probability of finding the ion in the state |6) is given by (see [Bastin et al. 2006]), 



|n = l([|6)(6|,i7o + if,.]) + (||&)(6|) 
in 



(81) 



The state of the ion at time t = to is givenby pni(^o) |+,/,) where p,„ is density matrix 
for the vibrational levels. Then the probability of finding the ion in state |&) at time t = io is 
given by, 

= Tr[p„(io)® (82) 



d 

d/' 



Therefore the generalized quadrature can be written in dimensionless time r = rjClt as, 
d 



(83) 

=0 

When choosing (j> — — 7r/2 and = the position and momentum operators are measured, 
respectively. Notice that only two measurements of the probability of finding the ion in state 
1 6) are needed to be taken to determine the position or momentum of the ion. 

In order to perform the simulation of the Zitterbewegung the ion can be cooled to the 
vibrational ground level, although it is not required, since it is also possible to measure the 
Zitterbewegung of a Doppler cooled ion. Then the internal state of the ion is set into \a). 
The next step is to excite the ion to a given motional state, i.e. a Fock state [Meekhof 
et al. 1996a, Meekhof et al. 1996/?], a coherent state [CaiTuthers & Nieto 1965], a thermal 
state [Stenholm 1986], a squeezed state [Heinzen & Wineland 1990, Cirac et al. 1993], or 
an arbitrary state. The relativistic-like dynamic ruled by the interaction Hjj (see equation 
Its! ) is generated by turning on the fields resonant with the blue and red sidebands at time 
t — 0. The time evolution of the ion ends at t — to when the interaction generated by the 
Hamiltonian Hu is turned off. Simultaneously, the internal state of the ion is driven to 1+^) 
by the application of an external field, then a JC coupling (see Hr in equation iTOli is applied 
to map the ion motional state into its internal state. At time t = to + t with r ^ 2tt/u!zb 
the population in the level |6) is measured by detecting resonance fluorescence. The previous 
sequence is rejjeated with the same parameters for a given number of repetitions to find the 
probability *(to + t). Then, the probability * at t = <o is also measured (without 
the need of applying the JC coupling). From the two last probabilities the expectation of the 
generalized quadrature can be deduced. The whole sequence is then repeated for different t 
to extract the trajectory of the ion. 
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6.4.3. Experiment of the simulation of the Klein 's paradox The Klein paradox [Klein 1929] 
is a relativistic effect described by the Dirac equation. It states that an electron can be 
transmitted with almost no reflection through a potential step with a height V bigger than 
twice the electron rest energy mc^. 

In the Dirac formalism the particle kinetic energy E takes into account the rest mass and 
therefore the electron momentum is given by 

p = -J{E-VY-m^c^. (84) 
c 

When — y| < m(? the momentum is imaginary and the wave function exponentially 
decays. However, when V > E + mc^, the momentum is real and one obtains a transmitted 
wave. This solution lies in the negative kinetic energy range E < —mc^ + V and therefore 
is associated with a electron-positron pair production. The positron sees an inverted potential 
step and therefore its motion is not impeded (see figure [T6]l. 




V+mc^ - 
V- 
V—mc^ - 




Figure 16. Potential Step (thick line). The white areas represent the forbidden regions where 
the solutions of the Dirac equation are exponentially decaying. The light and dark gray areas 
correspond to the positive and negative energy regions respectively. An electron with energy 
mc^ < E < V — mc^ is found initially in the positive energy region. For t > to, after the 
potential step is risen, the energy lies in the negative energy region allowing the creation of a 
hole in the Dirac sea. Picture adapted from [Schwabl 1997]. 

The transmission through a potential step can be reproduced by the Hamiltonian, 

i7« =i7«+y(|a)(a| + |6)(5|). (85) 

The last term is the potential step V{\a) {a\ + \b) {b\) that is set at the time t = tg. It can be 
implemented for example by the use of a detuned laser that produces the same stark shift in 
both levels. 

The experimental protocol is summarized in the following steps. First, the ions are 
cooled to the ground state \a). Then the JC and AJC couplings, described in equation [781 
are appHed at f = by turning on the red and blue sideband excitation fields. Next, the 
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potential step is turned rapidly on at time t = Iq. Last, the population in the state \h) is 
measured by fluorescence. The probability of finding the ion in the state | h) is obtained after 
several repetitions and is equivalent of the transmission probability through the step potential. 
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